Boost-invariant Hamiltonian approach to heavy quarkonia 
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Light-front Hamiltonian formulation of QCD with only one flavor of quarks is used in its simplest 
approximate version to calculate masses and boost-invariant wave functions of cc or bb mesons. The 
quark-antiquark Hamiltonian is obtained in the lowest (second) order of a weak-coupling expansion 
scheme for Hamiltonians of effective particles in quantum field theory. The derivation involves a 
heuristic ansatz for a gluon mass-gap that is meant to account for non-Abelian color dynamics of 
virtual effective particles in the Fock components with one or more effective gluons and can be 
improved order-by-order in the expansion. Fortunately, the resulting quark-antiquark Hamiltonian 
does not depend on any details of the ansatz within a large class of possibilities. It is shown that 
in the Hamiltonian approach in its simplest version the strong coupling constant a and quark mass 
m (for suitable values of the renormalization group parameter A that is used in the calculation), 
can be adjusted so that a) masses of 12 lightest well-established bb mesons are reproduced with 
accuracy better than 0.5% for all of them, which means 50 MeV in a few worst cases and on the 
order of 10 MeV in other cases, or b) masses of 11 lightest cc mesons are reproduced with accuracy 
better than 3% for all of them, which means better than 100 MeV in a few worst cases and on the 
order of 10 MeV in the other cases, while the parameters a and m are near the values expected in 
the cases a) and b) by analogy with other approaches. A 4th-order study in the same Hamiltonian 
scheme will be required to explicitly include renormalization group running of the parameters a and 
m from the scale set by masses of bosons W and Z down to the values of A that are suitable in 
the bound-state calculations. In principle, one can use the Hamiltonian approach to describe the 
structure, decay, production, and scattering of heavy quarkonia in all kinds of motion, including 
velocities arbitrarily close to the speed of light. This work is devoted exclusively to a pilot study of 
masses of the quarkonia in the simplest version of the approach. 

PACS numbers: 12.38.Aw, ll.15.Tk, 12.38.Lg 



I. INTRODUCTION 

In the Hamiltonian approach to QCD that is employed 
here, the calculation of masses of heavy quarkonia does 
not involve the well-known notions of scattering states 
for quarks and gluons, Feynman diagrams, path integral, 
euclidicity postulate, lattice gauge theory, or non-trivial 
vacuum expectation values. Instead, a renormalization 
group procedure for effective particles in quantum field 
theory (RGPEP, see below) is applied to quarks and glu- 
ons in canonical light-front (LF) QCD and leads to an 
effective boost-invariant Hamiltonian whose eigenvalue 
problem is expected to provide a first approximation to 
the true dynamics of the theory. The eigenvalues of the 
Hamiltonian are equal to the masses of the quarkonia 
(actually, squares of the masses) instead of their energies 
in a specific frame of reference. The approximate dy- 
namical picture studied here can be valid only for a set 
of states near the low end of the mass spectrum. 

The LF formalism is developed in a Fock space. Cre- 
ation and annihilation operators for effective quarks and 
gluons are calculable in a perturbative expansion using 
RGPEP (a brief review of the method is provided in the 
next section) and the basis states in the Fock space are 
formed by acting with the calculable creation operators 
on the state of vacuum. The effective Hamiltonian eigen- 
value problem exists in the momentum representation 
and would not be local if one attempted to write it in 
a position space. Nevertheless, the simplest approximate 



momentum-space dynamics respects not only the boost 
symmetry but also the rotational symmetry. Therefore, 
one can describe the structure of lowest-mass eigenstates 
corresponding to well-known hh and cc mesons using the 
spectroscopic scheme that is quite analogous to the one 
adopted in non-relativistic quantum mechanics with po- 
tential forces that respect rotational symmetry in a me- 
son center-of-mass frame of reference. Thus, the LF ap- 
proach has a potential to be helpful in solving conceptual 

Soblems with Poincare symmetry in quantum theories 
111 and the results described here can be considered 
an indicator of existence of a reasonable candidate for a 
new expansion method for solving theories as complex as 
QCD However, the approach is still in its infancy. 
The crude, heuristic study described here is merely a 
small step that needs to be taken on the way to find out 
if the RGPEP can work for heavy quarkonia in LFQCD 
as outlined in |j|. 
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Since the beginning of the theory of cc system H 
the de velopm ent of potential models 
sum rules ,2Qj , 

and studies based on effective Lagrangians |^ IH |23 , 
the latter method known also to work well in Hamil- 
tonian approach to QED p^. l25l |. current understand- 
ing of heavy quarkonia, especially in lattice approach 
in m m , is one of the best examples of g reat progress 
achieved in theory of strong interactions [29j . In the wake 
of the development, conceptual problems with a relativis- 
tic description of hadrons in the Minkowski space and 
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questions concerning the structure of the vacuum state 
continue to bother theorists 0,|33|- The problem is how 
to derive a quantum Hamiltonian for quarks and ghi- 
ons that precisely describes hadronic states in agreement 
with special relativity. 

One can write a formal expression for a Hamiltonian of 
quarks and gluons using various forms of dynamics 0. 
In the standard form, one describes the evolution of a 
system of particles using the time parameter t, which is 
measured along a time-like direction in space-time. But 
the formal expression must be regulated and renormal- 
ized, and one has to explain how to define the ground 
state of the theory that supports the tower of excita- 
tions that represent all kinds of relativistically moving 
and interacting hadrons. None of these tasks has been 
completed yet with the desired clarity and precision. In 
the LF form of dynamics, one uses a "time" parameter 
that is conventionally denoted by x"*" = t + z. Since the 
LF hyperplane a;"*" = is preserved by boosts along z- 
axis, and also by two other boost-like transformations, 
the LF form of Hamiltonian dynamics is invariant under 
the three special boosts. The three boosts are prerequi- 
site to the construction of hadronic states in all kinds of 
motion. The boost symmetry rises hope that the LF 
Hamiltonian dynamics can help in finding a universal 
theoretical description of hadrons in all frames of refer- 
ence, including both their center-of-mass frame (CMF), 
in which the constituent quark model is developed |3ll |. 
and the infinite momentum frame (IMF), in which the 
concept of partons is developed j32]. The exact descrip- 
tion of a hadron in motion is also essential in exclusive 
(or semi-exclusive) processes 33]. Regarding the vac- 
uum, the ground state problem in QCD does not appear 
in the LF form of dynamics in the way known from the 
standard approach |3(tI l33| and the LF Hamiltonian dy- 
namics of quarks and gluons is of great interest to many 
researchers [s^, as a serious alternative to the standard 
form. However, LFQCD challenges theoreticians with 
basic questions concerning quantum mechanics of parti- 
cles and fields and relativity. 

In search for understanding of the quark-gluon dynam- 
ics, it is natural to consider quantum chromodynamics of 
only heavy quarks coupled to gluons. A quark is consid- 
ered heavy when the phenomenological mass parameter 
associated with the quark, m, is much larger than Aqcd, 
the latter being defined in the RGPEP procedure that 
one can use to evaluate effective LF Hamiltonians '36^. 
The restriction to only heavy quarks creates a situation 
in which the renormalization group parameter, denoted 
by A, can be simultaneously much smaller than m and 
much larger than Aqcd- In such circumstance, the ef- 
fective color coupling constant at scale A, g\, can be for- 
mally considered small in comparison to 1 even when A is 
much smaller than m. The smallness of the coupling con- 
stant implies that the LF Hamiltonian of QCD expressed 
in terms of the effective quark degrees of freedom corre- 
sponding to A ^ m, denoted by H\, can be evaluated 
in RGPEP using perturbation theory. But the price to 



pay for this simplification is high because the quenched 
heavy-flavor version of the theory must be incomplete; it 
misses dynamical interplay among different flavors and 
entirely ignores effects due to light quarks. On the other 
hand, the price is worth paying because the single heavy 
flavor theory quickly renders a simple dynamical picture 
that may be helpful in learning more about LF QCD. 
Namely, it is sufficient to consider H\ obtained in just 
2nd-order perturbation theory and augment it with an 
ansatz for the mass gap for effective transverse gluons and 
these two steps already render a boost-invariant eigen- 
value equation in the effective quark- antiquark Fock sec- 
tor which has a well-defined and phenomenologically at- 
tractive structure 3 : the Coulomb potential with Breit- 
Fermi corrections is supplemented with a harmonic oscil- 
lator term with a frequency lu that is explicitly related 
to the values of a and m at the scales A at which this 
structure may be valid, which turns out to be the scale 
corresponding to the distances between quarks not much 
larger than the size of the lowest-mass mesons. At larger 
distances, additional gluons may be created and the po- 
tential may become linear, as will be discussed later. 

The resulting dynamical picture for the low-mass 
states is not sensitive to the details of the mass gap ansatz 
for effective gluons that was used to finesse the picture. 
Therefore, one needs to complete the 4th-order calcula- 
tion of H\ using RGPEP in order to begin a study of the 
true mass gap for gluons that may undergird the finessed 
effective picture. A major problem in the 4th-order cal- 
culation is to include light quarks. Even before the inclu- 
sion of light quarks, the first glimpse of the magnitude of 
4th-order effects could be obtained in one-flavor QCD by 
comparison of the 2nd-order picture with the 4th-order 
one. However, the 4th-order calculation requires detailed 
understanding of many complex terms in the LF Hamil- 
tonian at once even in one-flavor QCD. It thus becomes 
desirable to verify if the already known 2nd-order picture 
can provide a reasonably accurate description of data and 
become a candidate for the first approximation that can 
organize the 4th-order studies. 

Fortunately, in the limit of m large in comparison to A 
that itself is still much larger than Aqcd, one can con- 
siderably simplify the calculation of quark Hamiltonians 
using RGPEP. One can compute a 2nd-order Hamilto- 
nian that acts only in the effective heavy quark-antiquark 
sector, Hqq;^, and evaluate interactions in HqQy^ that 
depend on the quark spin. The result is that the lead- 
ing spin-dependent terms that one obtains from the 2nd- 
order RGPEP, in addition to the Coulomb potential and 
the harmonic oscillator force, automatically respect ro- 
tational symmetry. This leading simplest version of the 
Hamiltonian approach can be used to calculate the spec- 
trum of masses of cc or hh mesons. Comparison with 
experimental data is surprisingly optimistic. Not only 
the masses that one hopes to be reproducible relatively 
well turn out to be reproducible very well, but also those 
masses that one expects to be reproducible rather poorly 
turn out to be reproducible with considerably less accu- 
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racy than those that are reproduced well. Theoretical 
results seem to match data in a pattern that agrees with 
expectations concerning accuracy of the simplest version 
of the approach. 

Section briefly describes the derivation of the effec- 
tive Hamiltonians studied here. The leading approxima- 
tion to the heavy quark eigenvalue problem is explained 
in Section ITTll Section |Wl discusses spin effects and lists 
eigenvalue equations that describe mesons with spin 0, 1, 
and 2, with orbital angular momenta equal 0, 1, 2, or 3. 
Masses and wave functions of mesons that are obtained 
by numerical solution of these equations are discussed in 
Section Conclusions are summarized in Section IVII 
Key details of the calculations are relegated to Appen- 
dices. 



II. 



BOOST-INVARIANT HAMILTONIANS 



flow of the coupling constant. This effect becomes ex- 
plicit first in 3rd-order calculations [37l| . 
At a;+ = 0. the fermion field 
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and the gluon field 



are quantized by imposing commutation relations 



(7) 



(8) 



(9) 
(10) 



This section explains how the boost-invariant Hamil- 
tonians for heavy quarkonia that are used in the next 
sections to evaluate masses of cc and bb mesons are de- 
rived in QCD. One uses the gluon mass-gap ansatz in the 
intermediate step of evaluating Hqq^ from the 2nd-order 
Hx i. 



The measure of integration over momenta, [k], is 9{k^) 
dk^ dp k^ /(167r'^fc+) and the LF three-momentum 5- 
function is S'^ik^k') ^ 5{k+ -k'+) 5{k^-k'^) S{k^-k'^). 
Spins are denoted by a and colors by c. Further details 
concerning our notation can be found in the Appendix 
(see also Ref. |3|)- 



A. Canonical LF QCD 

One begins from the standard Lagrangian of color 
gauge theory with one flavor of quarks of mass m, 

£ = ^(z^ - m)^ - ^F^'-'F;, . (1) 

The corresponding generator of evolution in in the 
gauge — takes the form 

-ffcan = F[.^,2 + Ha2 + + Hj^i + H^Aip 

where each of the terms is an integral of a corresponding 
Hamiltonian density Ti. over the LF hyperplane, Hi ~ 
J dx^ (Px^T-Ci. Four terms that explicitly enter in the 
derivation of the approximate 2nd-order boost-invariant 
effective theory for heavy quarkonia, are 

1 - -9-"-^ -I- TO^ 

= 2^1^ Jq+ ^' 

Ha^ - ^^A^id^fA^, (4) 
W^A^ = g'ip^'ip, (5) 

H(^^)2 = ^g^^-f+t-^'jr^i^l^t-y^- (6) 

Other terms in the canonical Hamiltonian are also impor- 
tant. For example, the three-gluon coupling term plays 
an implicit role as a seed of the renormalization group 



B. Regularization 

The canonical Hamiltonian is divergent and the RG- 
PEP begins with regularization of the ultra-violet and 
small-A:"'' divergences. The ultraviolet divergences result 
from integrating over large transverse momenta of quanta 
in the intermediate states when one attempts to evaluate 
powers of the Hamiltonian. In 3+1 dimensional theory, 
the quadratically and logarithmically diverging trans- 
verse integrals result from momentum dependent spin 
factors for fermions and vector bosons. The small-A:"'' 
divergences arise due to gauge couplings of gluons. Note 
that in the yl+ = gauge only A^ and ijj^ — ^■f'^'f^ip 
are dynamical variables. In particular, A^ depends on 
A-^ and tp'^. As a consequence, interaction terms in the 
Hamiltonian can be written using the polarization vec- 
tor e** for gluons whose only two transverse components 
are independent degrees of freedom. A gluon with mo- 
mentum k~^ and k-^ has e~ — 2 e^k-^ /k^ . The fc+ in 
denominator in e~ is a source of small-A:^ singularities 
in LF QCD. There exist similar small-fc"*" singularities 
in the instantaneous interactions along z-axis on the LF 
hyperplane, especially where 1/9"*"^ appears, which hap- 
pens similarly to how the inverse of a three-dimensional 
Laplacian appears in the familiar Coulomb potential in 
the standard approach. The small-A:+ singularities also 
occur in a 1-1-1 dimensional theory |3^ . 

In the boost-invariant formulation of LF QCD, one 
does not regulate the theory by limiting the momentum 
components /c"*" and fc"*" of every individual quark and 
gluon separately. Instead, only the relative momenta 
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of particles in the interaction terms are limited. The 
regularization is accomplished by insertion of regulating 
factors r in vertices, see 0. The overarching rule for 
construction of the regulating factors r is that they must 
respect 7 kincmatical symmetries of the LF scheme. This 
rule is analogous to the requirement of rotational symme- 
try in the standard approach. There is a factor r = r/^rg 
for every bare particle in every vertex, ta limits the range 
of the relative transverse momentum of an interacting 
bare particle with respect to other particles participating 
in the interaction. The limit is set by the parameter A. If 
a particle of momentum k carries a fraction x of the total 
p"*" of particles in an interaction term, = xp'^ , then its 
relative transverse momentum with respect to the parti- 
cles interacting in this term is defined as k-^ = — xp^, 
where p-^ is the sum of transverse momenta of all par- 
ticles annihilated or (alternatively) created in the term. 
One uses ta = exp (— ^/A^), and the ultraviolet reg- 
ularization parameter A is sent to infinity in conimrison 
to all physical momentum scales, cf. [39ll40ll4ll [4^ . 

In the case of the small-fc"*" singularities, the regulat- 
ing factors rs limit the ratio x = k~^ /p'^ by the positive 
arbitrarily small dimcnsionless parameter S. All that is 
required of the factors rs is that they vanish as x^ when 
a; ^ 0. This condition is sufficient in globally colorless 
states. Linear divergences at small a; cancel out and one 
only needs to take care of the logarithmic divergences in 
integrals of the type J dx/x. The small- a: divergences in 
the gauge boson dynamics occur in both ultraviolet and 
infrared regimes. Massless particles can simultaneously 
have small x and small and their virtuality in the 
small- a; region can be large or small depending on the 
ratio of \k'^\ to ^Jx. Small x implies large virtuality only 
for particles with non-zero mass or fixed k^. 

Once the canonical Hamiltonian is regulated, one needs 
to introduce counterterms that restore the physics that 
existed outside the cutoff range. Thus, the counterterms 
remove effects of the regularization. For example, one 
inserts mass and vertex counterterms and they remove 
dependence on the artificial ultraviolet regularization fac- 
tors r A . The resulting regulated Hamiltonian with coun- 
terterms. 



H=[H,an+HcTUg. (11) 

provides the starting point for further steps. The further 
steps are also helpful in establishing the structure of the 
required counterterms. Note that the dynamics of color- 
singlet states of finite size should not be sensitive to the 
small- a; regularization. Namely, the singular limit a; — > 
concerns gluons with long wavelengths in the direction of 
x~ . But the strength of the coupling of such gluons to a 
finite-size color-neutral pair of quarks should disappear 
when the wavelength becomes infinitely larger than the 
distance between the quarks. 



C. Effective particles 

The initial Hamiltonian H of Eq. Hll(l is expressed 
in terms of the creation and annihilation operators de- 
fined by the Fourier components of local fields in Eqs. 
^ and lO. The same H can be expressed in terms of 
creation and annihilation operators for effective quarks 
and gluons that correspond to a renormalization group 
scale A in RGPEP. The procedure is constructed in such 
a way that the operators return to the canonical oper- 
ators when A tends to infinity, but when A is near the 
energy-scale of the binding mechanism, on the order of 
masses of hadrons, the operators create or annihilate ef- 
fective particles that are expected to correspond to the 
constituent quarks and gluons. The quantum numbers of 
the constituents are the same as in the local theory and 
one assumes that the corresponding creation and annihi- 
lation operators are related by a unitary transformation 

qx = UxqcanUl, (12) 

where the same letter q is used for both creation and 
annihilation operators. The next step is to express H in 
terms of q\ instead of qcam 

Hxiqx) = [Hcan + Hcrleg (Qcan) ■ (13) 

The Hamiltonian H remains the same but the coefficients 
in the expansion in powers of q\ are new. They include 
potentials whose structure can be calculated order-by- 
order in perturbation theory using the RGPEP. The key 
feature is that Hx has the structure 

Hx = fxGx , (14) 

where fx denotes form factors of width A and Gx repre- 
sents interaction vertices that can be calculated for any 
assumed shape of fx ■ The shape we use here is best de- 
scribed using the example of a term in which an effective 
quark emits an effective gluon: 

Gx = J[123]Gxil,2,3)al,bl^bx3, (15) 

fxGx = J[mfx{123)Gx{l,2,3)al,bl^bx3Xl6) 
/a(123) = exp [^{Ml^ - M^^X^] . (17) 

The invariant masses are defined by the formulas M12 — 
(/ci -f fc2)^, = /cf, using masses in the Lagrangian 
with 5 = to evaluate minus components of the four- 
momenta; m for quarks and for gluons. 

The operator Gx is defined by the coefficients of its 
expansion in a series of powers of operators qx ■ One can 
also define Qx , which is a series with the same coefficients 
but qx replaced by qcan- Then, one can use the constant 
operator basis qcan when solving differential equations of 
RGPEP for the coefficients, see js^. Qx is split into two 
parts: Qq and Qi = Q — Qq. Qo is the part that does 
not depend on the coupling constant g. The RGPEP 
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differential equation for Qj (prime denotes differentiation 
with respect to A) is: 

G'j - [fGi,m~m]'}gj , (18) 

where the curly bracket with subscript Qq is introduced 
to indicate the operator T that solves equation [T, Qq] = 
[(1 — f)Gi]' ■ The initial condition is that Qoc = H, and 
the solution is 

Gx = H + f ds [fsGls^m- fs)Gls\'}g,] ■ (19) 

This solution is evaluated order-by-order in powers of 
the coupling constant gx that appears in the vertices of 
G\ [s^]. The perturbative expansion is legitimate be- 
cause one never encounters small energy denominators. 
This feature is secured by the structure of the RGPEP 
equations and the shape of the form factor / (see the 
original literature). G\ is obtained from Qx by replacing 
Qcan by gA- 

Solving Eq. H19|l up to terms of order (there is 
no difference between the expansions in powers of the 
bare coupling constant g and the running coupling con- 
stant gx in 2nd-order terms, but one should think about 
the expansion in powers of gx, see below and next sec- 
tions), which includes finding the mass counterterms in 
the initial condition at A = oo, one obtains Hx that can 
change the number of effective particles by 0, 1, or 2. 
For 2nd-order evaluation of the effective Hamiltonian in 
the quark-antiquark sector, Hqqx, one only needs the 
following terms Q 

Hx = TqX + TgX + fx \Yqg\ + Vg^A + Zqqx] ■ (20) 

Tqx and Tgx denote the kinetic energy operators for 
quarks and gluons, respectively. Yqgx denotes the term 
which causes that effective quarks emit or absorb effec- 
tive gluons (the letter Y is chosen in the notation because 
its shape resembles an act of one particle splitting into 
two, or two particles forming one). fxVqqX is an inter- 
action between quarks due to exchange of gluons with 
virtuality greater than A. fxZqqx is the instantaneous in- 
teraction between effective quarks that originates in the 
instantaneous interaction in the canonical LF Hamilto- 
nian. Details of these terms are listed in Appendix IXI 

D. Derivation of Hqqx 

This section explains how one obtains the effective 
Hamiltonian Hqqx for a heavy quarkonium starting from 
the eigenvalue problem for the effective Hamiltonian Hx 
that reads 

Hx\P)^E\P). (21) 

|P) denotes an eigenstate of the operators and P^ 
with eigenvalues P+ and P-*- (see 0| for an example 
of RGPEP construction of the whole Poincare algebra). 



The eigenvalue has the form E = (M^ -I- P-'-^)/P+ and 
one obtains an eigenvalue equation for by multiplying 
Eq. (|^ by P+ and subtracting P^^. |P) is expanded 
in the effective particle basis as 

\P) = \QxQx) + \QxQxgx)+ ■■■ ■ (22) 

For A much smaller than m this expansion is dominated 
by its components with only two heavy quarks because 
the vertex form factors fx in Hx eliminate the probability 
of creating components with invariant masses that differ 
from 2m by much more than A. One may also expect 
that gluons develop a mass gap in QCD and the com- 
ponents with many gluons are also suppressed. If one 
neglects sectors with effective gluons entirely, the eigen- 
value problem is reduced to 

[TqX + fx {Vqqx + Zqqx)] \Q xQ x) = E\Q xQ x) • (23) 

But Eq. implies that the \QxQx9x) component sat- 
isfies equation 

[TqX + TgX + Vqqgx - E] [QxQxQx) = -Yx[QxQx) , (24) 

and can contribute to the dynamics in the sector [QxQx) 
in order 5^, or a a — 3a /(^^r), because Yx is of order ^a- 
VqqgX denotes potentials in the three-body sector, includ- 
ing non-Abelian potentials that act between the effective 
gluon and quarks. Additional interactions with sectors 
that contain four or more effective particles are not in- 
dicated. The additional interactions and Vqqgx are ex- 
pected to cause a shift in the gluon energy and make the 
eigenvalue equation differ from a similar one for positro- 
nium. In positronium, a state with two or more photons 
could have the same energy as the state with one photon. 
In QCD, there exist potential terms that act between glu- 
ons and quarks and among gluons themselves that have 
no counterpart in QED. It is very unlikely that there does 
not exist some shift in gluon energy that is absent in the 
case of photons. One can employ an ansatz for the effec- 
tive gluon mass in the three-body sector to study possible 
consequences of such shift . 

The point is that one can study the dynamics of Hx 
order by order in ^a using a scheme of successive ap- 
proximations that include an ansatz for effects that are 
extremely small for an infinitesimal ^a but need to be 
included to come close to a true solution that is obtained 
only when the coupling constant takes values compara- 
ble with 1. In each successive order one can replace the 
ansatz terms introduced in a lower order by a true inter- 
action of that lower order but with the coupling constant 
in them extrapolated to the large physical value 30] . The 
task of finding the initial ansatz terms that come close 
to the actual dynamics with large relativistic coupling 
constant may in principle require a lot of research to 
complete. Fortunately, the boost-invariant effective par- 
ticle approach has a useful feature: a lowest-order ansatz 
that is defined using an effective mass-like term for con- 
stituent gluons allows one to easily calculate and elim- 
inate the gluon component and the resulting dynamics 
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of constituent quarks comes out rotationally symmetric 
and independent of the details of the gluon mass ansatz 
one makes provided only that the ansatz satisfies some 
general conditions Q. In this first approximation, all 
interactions in the sector \Q\Q\gx) arc replaced by /i^, 
which is a function of the relative motion of the three con- 
stituents, /i'^ must vanish when the gluon a; — > 0. Since 
the mass ansatz for /i^ is supposed to model the domi- 
nant effect of all the interactions within the three-body 
sector and with sectors of larger numbers of constituent 
particles when the coupling constant takes a realistically 
large value, the first term in the ansatz can be considered 
to be on the order of 1 in comparison to the terms that 
depend on the infinitesimal coupling constant used in the 
formal expansion in RGPEP. The whole eigenvalue prob- 
lem for H\ with the ansatz is now reduced to only two 
coupled equations (A is omitted). 



Y\QQg) 



{T, + fg)\QQg)+Y\QQ) 



E\QQg)l25) 
E\QQ). (26) 



The operator Tg is marked with the tilde in order to 
indicate that the effective gluon mass in Eq. (jA3|) is 
replaced by the ansatz mass /i^ in the 3-body sector. 

The Hamiltonian Hqq that acts only in the \QQ) sec- 
tor can now be evaluated as a power series in g using 
an operator usually denoted by R In the simplest 

version, R expresses the 3-body component through the 
2-bodyone, \QQg) = R\QQ). Note that the sector \QQg) 
is separated from the sector \QQ) by a gap in invariant 
mass. The 2nd-order result is the Hamiltonian whose 
matrix elements are ^ 



{1?.\Hqq\2A) = (13| [T, + / {V,, + Z,,)] |24) 
1/2 1/2 



(13|/n, 



£^24 — Tq — Tg 



E 



13 



T„-r„ 



where with i equal 1 or 2 and j equal 3 or 4 are 
eigenstates of the operator Tq in the \QQ) sector of the 
Fock space, and Eij are the corresponding eigenvalues. 
The labeling of states is illustrated in Fig. [3 The basis 
states are defined as 



= 44, lo), 



(27) 



where h\ and d\ are creation operators for effective 
quarks and antiquarks corresponding to the RGPEP 
width parameter A. The corresponding eigenvalue is 
E,j = {Ml+P^^)/P+, wheveMl = i^. + k^f and the 
minus components of the four-momenta are evaluated as 
for free particles of mass to. 

In summary, the procedure used here p| replaces the 
eigenvalue problem for H\ by an eigenvalue problem with 
an ansatz (dots denote operators that couple states with 
more effective particles than three) 



H3 Y 
H2 



T3 + ^i' Y 
Ft T2 + V2 



,(28) 



and then the operation R is used to derive the effective 
quark Hamiltonian (in a simplified symbolic notation) 

1 



H, 



QQA 



= T2X + V2X + Y, 



A 7^3 



;Yx. 



(29) 



The procedure should not be confused with a conven- 
tional Tamm-Dancoff approach to quantum field theory. 
The quantum particle degrees of freedom that are ob- 
tained from RGPEP are not the bare quanta of local 
canonical theory, see |46(| , and the effective particles obey 
rules of the LF dynamics with a vacuum that is simple 
to work with. Moreover, the effective particles interact 
through terms like Yx that contain vertex form factors 
whose width has interpretation of the size of the effective 
particles in strong interactions (the particles cannot emit 
or absorb any quanta with greater invariant mass changes 
than A). At the same time, A also plays the role of the 
RGPEP parameter in the differential equations that con- 
trol the evolution of operators from the canonical ones at 
A = 00 to the effective ones that can be used in a rela- 
tivistic computation of bound states when A is lowered to 
the scale of the hadronic masses. The relativistic nature 
of the procedure is reflected by the possibility to con- 
struct all generators of the Poincare group at the scale A 
at which one wishes to solve the eigenvalue problem for 
Hx itself "il. 

Although the eigenvalue problem for heavy quarkonia 
in QCD with explicit inclusion of the sector with 3 ef- 
fective particles has not been studied in detail yet, it is 
important to state here that the ansatz scheme dictates 
in this case the replacement 



H4 Yi 
Y^ H3 

yI yt U2 



Y2 

Y 



Yi 

1 H: 



Y2 
Y 



(30) 



and subsequent application of R to the desired order. 
These steps appear to resemble the LF Tamm-Dancoff 
scheme with sector-dependent counterterms proposed by 
Perry, Harindranath, and Wilson '47'|. The concep- 
tual difference is that in the Perry-Harindranath- Wilson 
scheme the elimination of sectors occurs within a Hamil- 
tonian eigenvalue problem with large cutoffs and the 
ultraviolet renormalization issue is a part of the prob- 
lem, leading to Wilson's triangle of renormalization with 
a vast space of relativistic quantum operators. In the 
scheme used here the ultraviolet renormalization 
group procedure is completed long before one tackles the 
eigenvalue problem and introduces an ansatz at the scale 
A near the magnitude of invariant masses that character- 
ize observables. The principles of extracting a small and 
computer-soluble eigenvalue problem from an eigenvalue 
problem of infinite size, such that the results obtained 
from a small problem can represent solutions to the infi- 
nite problem, will be further explained below. 

Finally, it should be stressed that the mass ansatz has 
the structure 01 



Ma 



[I - al/a^^ /i' 



(31) 
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where as denotes the large, relativistic couphng constant 
that the effective couphng ax is supposed to reach as a 
result of extrapolation from the infinitesimal values used 
in the perturbative solution of the equations of RGPEP 
for Hx. This structure ensures that in the lowest-order 
approximate expressions obtained through expansion in 
powers of aA, only the first term counts, 



2 

l^ansatz 



(32) 



This term is order 1 . But when one increases the order of 
included terms and extrapolates to ax = cts , the ansatz is 
removed and one has a chance to recover a true solution of 
the initial theory with increasing accuracy. Comparisons 
of results obtained from expansions in successive orders, 
and use of better extrapolation techniques than through 
a plain power series, will hopefully indicate if the proce- 
dure can converge on a well-defined dynamical picture. 
For that purpose, one should compare the initial approx- 
imate theories with experiment and find out if the cou- 
pling constants required for explaining data can be small 
enough to pursue the chain of calculations based on the 
weak-coupling expansion. Results of this study suggest 
that at least for heavy quarkonia the required coupling 
constant in LF QCD is considerably smaller than 1, see 
below. 



E. Eigenvalue equation for Hqq^ 
The eigenvalue equation for Hqqx has the form 

HqQx\QxQ>.) = ^ \QxQx). (33) 

The eigenstates are written as (see Appendices H and |b1 
for details of the notation, subscript A is omitted) 



M~5P+^^Vj). (34) 
— 4m^ 



The eigenstate wave function can be written in the 
form that exhibits its covariance under 7 kinematical LF 
Poincare transformations, 

^ ^^'^ s,Sj{kij)uk^^SiVkj,Sj , (35) 

SiSj 

where ^ cMFij{kij) denotes the wave function that de- 
pends on the relative three-momentum of quarks in their 
CMF, assuming their masses are just m. The indices 
Si and Sj denote projections of spin on z-axis. Spinors 
Si E^nd Wfe s are obtained using LF boost matrix 
(A±=7V/2)' 



B(k, m) 



1 



[k+K+ + K-{m + k^a^)\ (37) 



acting on the spinors at rest, itg^s and tiQ^s, in the refer- 
ence frame in which the bound state calculation is carried 
out and where the four-momentum of the bound state has 
components P+, P^, and P" = {M"^ + p^'^) / P+ , M be- 
ing the eigenvalue that one wants to calculate. Spinors 
u^: ^ and v_^_ ^ are obtained by "boosting" spinors at 
rest in the CMF of the constituent fermions along their 
relative three-momentum. An additional spatial rotation 
is applied to spinors in the CMF before the latter boost 
is applied, in order to build a spin basis in which one ob- 
tains explicit rotational symmetry of spin-dependent in- 
teractions in the leading approximation. The additional 
rotation is the same as the well-known Melosh transfor- 
mation [isl I49I I . Details of our notation for momentum 
variables, spinors, and boost matrices, are explained in 
Appendices El and O 

Eq. H33(l implies the eigenvalue equation for the wave 
function, 



K 



13 



m 



Yl 



a;3 



M 



/f2 



*si.3(fcl3) 



167r3 



X2X4^ 



^«A(13,24)vI/,^,Jfc24), (38) 



in which the mass- like terms rnyi and result from 
the self-interaction of effective quarks through emission 
and re-absorption of effective gluons, and ua(13, 24) re- 
sults from the exchange of the effective gluon between 
the two quarks. The ultraviolet- finite part of the mass 
counterterm in the effective quark mass m\ is so adjusted 
(using a single quark eigenvalue problem) that at A = Ao 
one obtains B 





X1X3 


Xl 


la' f 


dx2(pK24 


167r3 J 


X2XA 



X3 



*siS3(fcl3) 



^i;o(13,24)*,,,,(fc24) =0, 



(39) 



where 



t;o(13,24) = -Ag,^j^,r,, + B^-f^, (40) 

and Snii, (5m|, A, P, and other symbols, are explained 
in detail in Appendix IbI 



III. 



LEADING APPROXIMATION 



The RGPEP result of Eq. ^ is further analyzed 
as a typical window Hamiltonian eigenvalue problem of 
the kind studied in detail in the case of a generic ma- 
trix model with asymptotic freedom and a bound state 
[sol Isil Is^ . The model is soluble exactly and provides 
a relatively well-understood pattern to follow in the case 
of QCD with one heavy flavor. Earlier LF studies, based 
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on coupling coherence [s^, [3 EE EE E3 1 did not have 
such a pattern to follow and did not use a boost-invariant 
concept of effective particles. They were carried out in a 
frame of reference nearly at rest with respect to the CMF 
of the quarkonium, and employed a logarithmically con- 
fining potential that was obtained in the quark-antiquark 
sector neglecting all other sectors. Those studies pio- 
neered an attack on the bound-state problem in LF QCD 
along the path discussed by Perry [5J, |53 , including el- 
ements of the method outlined in Wilson et al. |30| . 
such as the absolute cutoffs on momentum variables (es- 
pecially fc+) in the Fourier analysis of field variables in 
position space, or similarity RG procedure for Hamiltoni- 
ans. One of the key issues of the LF approach is how to 
obtain rotational symmetry and the initial studies had 
to struggle with the issue, in addition to the issue of 
construction of counterterms that restore boost symme- 
try violated by cutoffs on absolute momentum variables. 
The RGPEP procedure used here leads in its simplest 
version to a boost-invariant and rotationally symmetric 
spectrum of meson masses. 



A. Coupling constant in the window 

Construction of a window eigenvalue problem begins 
with a selection of a set of states of effective particles 
with kinetic energies (actually, free invariant masses) in 
a certain range that is also called a window, for brevity. 
The size of this range should be larger than the width Aq 
which appears in the form factors / in Hxg ■ 

The next step is to evaluate matrix elements of 
in the selected window of basis states. These matrix 
elements form a matrix W of the window Hamiltonian 
whose non-perturbative diagonalization is to produce the 
bound state of interest. To facilitate efficient diagonal- 
ization in a continuum theory, one can use a set of or- 
thonormal wave packets (such as the wave functions that 
solve a two-body bound-state eigenvalue problem with 
a harmonic oscillator potential) as a basis in which the 
matrix elements of the window are evaluated. 

Typically, if the energy range (the word "energy" 
should be replaced by the words "invariant mass" almost 
everywhere in this paper, but the reader is expected to 
be more familiar with the word energy than invariant 
mass in reference to the quantum dynamical concepts 
that count here and we use the word energy in order to 
avoid confusion due to the lack of familiarity with the LF 
form of quantum Hamiltonian dynamics) in the window 
Hamiltonian matrix W is sufficiently larger than Aq, the 
middle eigenvalues of the window are independent of the 
window boundaries and they match the eigenvalues of the 
full H\„. The latter eigenvalues are equal to the exact 
eigenvalues of the initial Hamiltonian H if the RGPEP 
procedure is carried out exactly. The additional virtue of 
lessons from Ref. [s^ . beyond showing that an asymp- 
totically free model can be solved using a window, is that 
one can also evaluate the matrix elements of W in pertur- 



bation theory, as if the coupling constant was extremely 
small. One sets the coupling constant to a realistically 
large value when one solves the non-perturbative eigen- 
value problem for W. 

The point is that a few low orders in RGPEP calcula- 
tion of W may lead to a good approximation (reaching 
better accuracy than 10% for the bound-state eigenval- 
ues already when W is calculated in second order) if one 
properly chooses Aq in order to work with a small number 
of basis states (small means small enough so that they 
can be handled using computers) and if one adjusts the 
coupling constant in the window to the chosen Aq. Aq 
should be near the scale of invariant masses that domi- 
nate in the binding mechanism. The coupling constant 
is defined through the value of a specific matrix element 
in the window. It is adjusted by comparison of the spec- 
trum of W with data (in [s^l, the role of data is played 
by the known exact spectrum). The main result of the 
matrix model (studied so far up to 6th order, or 5 loop 
integrals) is that when the coupling constant in the win- 
dow W is adjusted so that one middle eigenvalue of W 
matches the corresponding exact solution then also other 
middle eigenvalues of W approximate the corresponding 
exact solutions. 

In a theory as complex as QCD, the RGPEP calcu- 
lations of window Hamiltonians beyond 2nd order will 
require a lot of work. Completion of the 4th order calcu- 
lation is important for determination of the accuracy one 
can achieve using window Hamiltonians in QCD. Apart 
from the plain perturbative expansion, one may eventu- 
ally take advantage of the idea of coupling coherence EE] , 
reformulated for the use in RGPEP. However, when ad- 
ditional flavors of quarks are included and their masses 
are lowered toward small values known in the standard 
model, one may have to deal with an infrared limit cycle 
EEl (and universality that may apply in that case |59j . 
instead of the asymptotic freedom structure known in the 
ultraviolet). But in the case of heavy quarks, i.e., when 
the quark mass is formally very large in comparison to 
Aqc d , the value of the coupling constant required in the 
window may be small and no complications possible for 
light quarks are expected to occur. 



B. Heavy quark limit 

In a formal analysis of QCD with one heavy flavor, 
the quark mass m can be much larger than Aq and the 
latter much larger than Aqcd- In these circumstances, 
the perturbative coupling constant corresponding to Aq 
is small and the relative motion of quarks in the sector 
\QQ) is limited by the form factors / of width Aq. This 
means that the dominant relative momenta of quarks in 
mesons are small in comparison to m. (One should ob- 
serve that the process of extrapolation of the window 
to a large value of the coupling constant corresponds to 
the increase in the value of Aqcd', at a phenomenolog- 
ically useful value of ao = aAo, one may have to work 
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with Ao that is no longer very small in comparison to 
m). Formally, in the limit of large to, the dominant dy- 
namical effects in the window eigenvalue Eq. 1)39(1 with 
Xo ^ m can be analyzed using a non-relativistic approx- 
imation for the relative motion of quarks. The arbitrary 
motion of the meson as a whole can be exactly sepa- 
rated from the relative motion of the constituent quarks 
because the RGPEP is invariant under the LF boosts. 
Therefore, approaching the IMF, one will control the dif- 
ferences between absolute values of momenta of the quan- 
tum constituents that extend to much larger values than 
the mass of the entire meson. The differences correspond 
to some fixed range of Feynman x around 1/2. But every 
fixed value of x different from 1/2 implies that the corre- 
sponding quark momenta in the IMF differ by amounts 
infinitely larger than any fixed mass parameter. 

The non-relativistic approximation is formally vali- 
dated by the condition that Aq <C to, and that the in- 
teraction terms are growing not faster than a polynomial 
of kinetic energy (invariant mass) and cannot overcome 
the exponential falloff of the form factors / for changes 
of invariant masses that are larger than Xq. Thus, in or- 
der to use the non-relativistic expansion, one has to keep 
in place the exponential form factors / that provide the 
convergence - these form factors determine the size of the 
window in momentum space and one cannot expand them 
in a series. What can be expanded are the arguments of 
the exponential functions, the perturbative factors that 
appear in RGPEP in addition to the form factors, spin 
dependent factors in interaction vertices, and relativistic 
measure of integration over relative momenta of quarks. 
The accuracy of results of diagonalization of the window 
W\g will depend on the choice for the form factor func- 
tion / and, especially, on the optimization factors that 
are critical for the convergence of the perturbative eval- 
uation of H\g in 4th order |^ The optimization 
factors were considered in the case of heavy quarkonia 
elsewhere (Gfll • The study described here was carried out 
using / of the generic type exemplified in Eq. ((T7|) . Every 
form factor considered here is the same exponential func- 
tion of the square of a difference of squares of invariant 
masses of the effective particles in interaction. 

When Eq. (|39|) is written in the non-relativistic ap- 
proximation for the RGPEP, spin, and integration mea- 
sure factors, the limit of small coupling exhibits scaling 
property similar to the Schrodinger equation for positron- 
ium. The scaling in RGPEP is described at the beginning 
of Appendix^ The scaling implies that the quark eigen- 
value problem is dominated by the relative momenta on 
the order of strong Bohr momentum, 

7 4 TO 

kB = 3"oy, (41) 

which is the quark analog of a me/2 in positronium in 
QED. This scale emerges from the window eigenvalue 
condition independently of the value of Aq as long as 
Aq is sufficiently large in comparison to kg. In other 
words, the eigenvalues Af^ depend mainly on the value 



of and not on the value of Aq itself when Aq is much 
larger than ks- The questions of how large Aq must be 
in comparison to ks in order to obtain results that are 
sensitive to Aq practically only through the value of ao, 
or to what extent this lack of direct sensitivity to Aq 
itself is obtained for realistic values of ao and Ao, are 
not answered here. A study of such issues has been done 
before in a model based on Yukawa theory 46, 61] . 

In a formal analysis of the non-relativistic expansion 
for infinitesimal coupling constant ao, one can assume 

TTl 

Ao = A„y, (42) 

Am = f-aoj Xp. (43) 

Thus, Ao is much smaller than to/ 2 when Ap is on the or- 
der of 1. The two parameters Ap and e are useful in sep- 
arating different terms in the complex, spin-dependent 
interactions that otherwise do not occur clearly ordered 
in size. The particular choice of the power 0.5 -I- e follows 
from how the form factors / limit momentum transfers 
in vertices. The most interesting case is e close to (see 
below). At the same time, kg is much smaller than Aq 
as long as < e < 0.5. In the formal analysis, kg is con- 
sidered much larger than Aqqd when one evaluates the 
window Wx„ using RGPEP. But when one extrapolates 
the coupling constant in the window to realistic values, 
a realistically large value of Aqcd is introduced, instead 
of an infinitesimally small one. The ansatz for the gluon 
mass gap corresponds to the scale of the realistic Aqcd ■ 
Note that the formally introduced relationship between 
Ao and ao does not mean that one replaces the true RG- 
PEP dependence of ao on Ao by an artificial one. All 
that is done is to introduce a parameterization of an 
unknown infinitesimal value of ao at a single value of 
Ao; the parameters A^, Ap, and e, remain free to change 
while kg stays always much smaller than Ao. After the 
scaling picture is described using this parameterization 
and identifying terms that scale with different powers of 
ao, one can look for the values of to and Am for which 
the scaling picture extrapolated to large values of the 
coupling constant is useful phenomenologically. At that 
point one identifies the realistic values of ao and Ao. All 
one obtains this way is an approximate picture for heavy 
quarkonium dynamics that can serve as a starting point 
for a systematic calculation of corrections using RGPEP. 

C. Scaling of different factors 

The scaling expansion of RGPEP factors A and B that 
occur in the potential fo(13, 24) in Eq. H39|l . in terms of 
powers of ao (we shall omit the subscript below when- 
ever it is irrelevant to the context) is given in Eqs. (|B43() 
and (|B44|) . The first terms in the scaling expansion that 
provide the leading approximation are 

A^~f^ + ff^^^, (44) 
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B 



// 



Qz 1^ + 



(45) 



The symbol q denotes the three-momentum transferred 
between quarks. The self-interaction terms Sm? contain 
the same function B in their integrands. 
The RGPEP factors // behave as 



// - exp 



/ mq^ 



(46) 



cancel out or get corrected when the window is calcu- 
lated in 4th-order RGPEP. Since the spin factors have 
distinct origin (quark current factors that are specific to 
QCD) than the generic RGPEP factors (the same in all 
field theories) and momentum space integration measure 
(universal in relativistic particle physics), the structure 
of interaction kernels that one encounters in the scaling 
equation (in addition to the spin-independent harmonic 
oscillator term) can be written as (the factor in front 
is not included) 



Using definitions of Eqs. (jB25|) and ljB26p for momenta 
and and (|^ for Aq, one obtains 



exp ■ 



8p2 



[(|a)'Ap] cos2 0. 



(47) 



where q^ was replaced by q cos 9 and 9 is the angle be- 
tween q and z-axis. For small a, the form factors are 
not zero only for small p, and, in fact, only for vanish- 
ingly small p when the angle 9 between p and z-axis ap- 
proaches Tr/2. Unless there exists a large contribution 
in the region of small p, especially near 9 ~ 7r/2, from 
other factors in the potential, the factor // is equivalent 
to zero when a is near in the scaling analysis. But B 
can be large for 9 7r/2 due to q^ in denominator. This 
singular behavior originates in the instantaneous LF po- 
tential due to gauge coupling between colored particles. 
One has to find the result that survives in the limit of 
small a in the presence of the singularity. The factor 
regulating the singularity at qz — is provided by the 
ansatz function /u^. When one combines the terms with 
// in the self-interactions Sm^ and in A and B in the 
gluon exchange potential, the net result is a spherically 
symmetric and spin-independent harmonic oscillator po- 
tential whose spring constant is no longer sensitive to the 
mass ansatz /i^ under quite general assumptions 4] . The 
oscillator frequency is 



4 a 
3 n 



A 



V1152 



1/4 



(48) 



and the corresponding spring constant, k — muj'^ /2, leads 
in the dimensionless Schrodinger equation in variables 
Pij, defined in Eq. (|B25II . to the oscillator term that 
scales like a^*^ and becomes independent of a when e — > 0, 
see Eqs. (|50f) to (|53|) below. The oscillator potential is 
independent of the quark spins. 

According to Eq. H40|) , all spin effects in the leading ap- 
proximation originate from the current factors that mul- 
tiply the term —f/q'^ in A in Eq. H44|l . The leading spin 
effects can be identified using the same scaling analysis. 
In the scaling analysis for infinitesimal a, spin-dependent 
terms are times smaller than the spin-independent 
terms. Therefore, one can also try to include correc- 
tions order that do not depend on spin. However, the 
spin-independent correction terms order a"^ that emerge 
from the scaling expansion based on only 2nd-order RG- 
PEP violate rotational symmetry and are expected to 



V 



{l + a'^S){l + a^M) , (49) 



where R refers to the RGPEP factors, S to spin, and M 
to measure. 

We drop the term R because it depends on the z-axis 
and can only be corrected in the 4th order calculation 
for RGPEP factors. In the 2nd-order RGPEP factors 
analyzed here, the correction R is given by the term c 
in Eq. ljB43p and when one averages this term over the 
direction of the z-axis, it vanishes. Nevertheless, one 
should remember that a genuine 4th-order calculation 
of H\ may produce corrections of the type R that will 
change the radial dependence of the potential from the 
Coulomb shape to a different one. 

In order to identify the leading spin effects, we combine 
the spin and measure corrections to the factor l + a'^{S + 
M) and write it shortly as 1 -I- BF, where BF stands for 
Breit-Fermi terms. The point is that one can observe 
cancellation between S and M factors and the remaining 
terms produce a rotationally symmetric spin-dependent 
terms after one introduces the additional turn in spin 
basis that is described in Appendix lO 



D. Structure of the eigenvalue problem 

Finally, using dimensionless variables pi^ and p24 (Fig. 
Elillustrates the labeling of the momentum variables) de- 
fined in Eqs. ljB25|l and ljB26|l . one arrives at the fol- 
lowing eigenvalue equation for the spin-dependent 2x2 
matrix wave function (j), defined in Eqs. (|D4p and ljD5|) : 



= [Pi3 - K \i3 - 0(Pi3) - 2 



(27r)3 



V 



47r , 
/ — (I 



BF), 



f = exp <^ - 



—a 
3 



l-2e 2 
Pl3 



H 2 ' 



A^ 



krt — 



VlT52^ 16 V3" 



'^(t>{P2i) , 

(50) 
(51) 

(52) 
(53) 



In the limit a — > in the above result, the eigenvalues x 
tend to —l/'n? with natural n (the Coulomb spectrum). 
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One obtains meson masses by evaluating the eigenvalues 
X for realistic values of a and using the formula 



IV. 



SPIN EFFECTS 



M = 2m\ll + x \ -a 



(54) 



The fact that a mass-gap ansatz for giuons leads to 
an oscillator-like interaction term, which respects rota- 
tional symmetry already in the 2nd-order analysis, in 
which fj,^ ^ 1, does not seem accidental. The result is 
almost independent of all details of the ansatz because 

is limited by the form factors // in the function B of 
Eq. (|45(l to so small values that the ratio n'^ /{q^ + A*^) 
is practically 1 for any reasonable ansatz. In addition, 
it seems not excluded that the same result comes out 
also as a part of the genuine 4th-order calculation. In 
the 4th-order calculation, the part of the ansatz for fi^ 
that is order 1 may cancel out in the window W with 
large a. Looking at Eqs. H29|l and 129|l . one can see 
that the ansatz term order 1 in the three-body sector 
is eliminated in 4th-order RGPEP calculation when one 
includes the three-body sector in the non-perturbative 
window dynamics. If instead one uses the perturbative 
operation R to further reduce the window to the two- 
body quark-antiquark sector only, the cancellation of the 
ansatz begins in 6th-order calculation. For large cou- 
plings, there can exist cancellations that cannot be easily 
foreseen. But even if the ansatz is not introduced at all, 
some shift of the three-body invariant mass, say /ij, will 
emerge from QCD interactions of formal order a in the 3- 
body sector and this is how the actual gap may show up 
for the first time. When one proceeds to the scaling anal- 
ysis of functions analogous to A and B in Eqs. H44|l and 
(ED), the new shift will be of order a if it is proportional 
to rn^, of order a'^/^ if it is proportional to mA, and of 
order if it is proportional to A^. But the momentum 
transfer squared, q^, is of order and it may continue 
to be formally much smaller than or comparable to /i^ 
in the scaling limit of small a. The form factor / limits 
|g I to values on the order of (|i|/2) A^ (4a/3)^'^ ks, where 
t — cos 9 and 9 is the angle between the momentum of 
the effective gluon and the z-axis. The terms that lead 
to the harmonic force originate from the singular behav- 
ior of q~^ ^ when t ^ 0. But one can still neglect 
q^ ~ t^a^+'*' in comparison to ^\ ^ a^'t^'^^" with n — 1 
and 3/2, and even for n = 2 the result of integration 
may be close to the one obtainable when q^ is neglected 
in comparison to /i^, cf. Q. Also, if the effective mass 
ansatz is just a first term in the expansion of the gluon 
gap in powers of the gluon momentum squared, which 
corresponds to the case with ~ q^t'' in the limit of 
small t 4] , the scaling applies in the same way and leads 
to the same oscillator result 0- So, if some mass gap 
shows up in order a in the 3-body sector, as one expects 
it to happen in QCD, the results of Eqs. and 
may still be valid. 



Spin effects are caused by the Breit-Fermi terms, BF 
in Eq. H51I) , that originate from the product of currents 
j 12 jiSfi ^0 in Eq. (|40|) . In the leading approximation, 
vq is displayed in Eq. (|Dip . The spinors Uk^.si and Vk-.s- 
in the currents originate from the canonical LF Hamilto- 
nian of QCD and they are related by boosts Cij described 
in Appendix [HI to the spinors ^ and ^ that are 
introduced in the definition of the CMF wave function 
^CMFijikij) in Eq. 

It is shown in Appendix IdI that the BF terms in Eq. 
(|ET)l are 



BF^^- 



Kpli + Piz) a - paba pa 
Spi^a p24<^ ba + 3b(j p24<^ Pi3<^ , (55) 



where (j) denotes the 2x2 matrix wave function cj) — 
a + ba and p is the difference between pi^ and p24- An 
alternative form of the same result, 



BFcj) 



(P24 



a' 
"3 

"9 

b (8p24 



Pis) a 



(4pi3P24 +P?3 +P24) ba (56) 
2pi3) - b (4pi3 + 2p24) P24<? , 



shows that the singlet wave function a and the triplet 
wave function b are not coupled and describe different 
eigenstates. 

The resulting eigenvalue equations for different mesons 
are listed in subsections below. These equations are 
boost-invariant and describe the relative motion of two 
heavy quarks no matter how fast the whole quarkonium 
is moving, which is also reflected in Eq. H54|l that differs 
from a non-relativistic expression for energy of a slowly 
moving object, E = M + Eb + P^/(2M), where the 
binding energy Eb is given by some Schrodinger equa- 
tion. In order to obtain a state of a moving quarkonium, 
one has to insert the wave function (j> into Eq. H34() . us- 
ing Eqs. 1(23), (ISnil, EH), and ^H^, all of which are 
relativistic. Note that the eigenvalue equations in the 
subsections below lead to non-local interactions at short 
distances if one introduces position variables canonically 
coupled with the relative three-momenta of quarks de- 
fined here in a boost invariant way. The reason for non- 
locality is that the interactions contain the form factors, 
/, that depend on the differences of invariant masses be- 
fore and after an interaction. One should remember that 
although the equations in the next subsection look decep- 
tively simple and similar to non-relativistic models, they 
describe wave functions that from the models' point of 
view formally correspond to different frames of reference 
for different values of the modulus of the relative three- 
momentum of quarks, see Appendices. 
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Note also that the interaction terms in the equations 
Usted below include potentials that in the absence of the 
form factors / would produce three-dimensional Dirac 
(5-functions in the position space representation. The 
(5- functions would lead to ultraviolet divergences. But 
the form factors / smear the (5-functions and render fi- 
nite results. Nevertheless, the terms with the smeared 
(5-functions are not weak when a is extrapolated to val- 
ues on the order of 1 and they contribute to significant 
spin effects for states that involve significant s-wave com- 
ponents. 



A. Mesons ric and rji, 

In this case, 6 = and the eigenvalue equation for 
4'iPij) = '^ij/Pij takes the form 

= [pf3-fc.A,,3-x]^- /|5iv^,(57) 



2p 



13 



(2^)3 P24 



where 



1 + Yb24+P?3) 



(58) 



The orbital angular momentum is zero and the inte- 
gration over angles (see Appendix Ie)) produces a one- 
dimensional integral equation 



= h. 



: 013 



2 f°° 

- / dp24 / P13P24 VV a24 , (59) 

I" Jo 



with 



1 + Y(P24+P?3) 



■^0 J 



where the function Jq is given in Appendix IeI 

hsosc — Pi3 ^p'5l3 2; 



(60) 



(61) 



is introduced as a generic notation for the s-wave har- 
monic oscillator terms in the mass eigenvalue equations 
for all mesons. 



B. Mesons ,//* and T 

In this case a = and the eigenvalue equation describes 
a function (p{pij) — bij<7, where 



D satisfy two coupled equations 

^.so.sn, 







^SOSC 1 

0, h 



sosc ~t- kp 



Pl3 



Sl3 
Di3 



dp24 f P13P24 







S24 






D24 



where 



Wss = Ja + Y [(^13 + P24) ^0 - 16/9] , 
Wsd = y [p?3 iJ2 - ^0) + 4/3] ^ , 



Wdd 



[pL {J2 - Jo) + 4/3] 



V2 
3 ' 



J2 + (J2 - Jo)/2 

^,2 



(63) 

(64) 

(65) 

(66) 
(67) 



+ Y {(P?3 +P24) [^0 - (J2 - Jo)/6] - 20/9} , 
and the functions Jq and J2 are given in Appendix IeI 

C. Mesons Xcq and Xbg 

Here a — and the eigenvalue equation for 4'{'Pij) = 
bij pij a /plj takes the form 

= ( hsosc + kp^\ 613 - - / dp2i f P13P24: W 624 : 



Pl3 



where (see Appendix lEjl 

^,2 



(68) 



W = Jl + y [P13P248J0 - (p?3 +P24) ^1] ■ (69) 



D. Mesons Xci or Xbi 

Here again a = and the radial eigenvalue equation 
for (pipij) — bij sx Pij <^/Pij takes the form 

/ 2 \ 2 /"^ 

hsosc + kp^ ] bi3 / dp24/Pl3P24 ^624 , 

V Pl3 / ^ Jo 



(70) 



where 



W = Ji + Y [2Pi3P24(Jo + J2) + {pIs+pIa) Ji] 



(71) 



"13 



srkl ^13 



1 



Pl3 V2 



o Pl3H3 \ -P13 
Pl3 ) Pl3 



(62) 



and s is a polarization vector of a massive meson of spin 
1 . The s-wave wave function S and d-wave wave function 



E. Mesons Xca or Xbs 



In this case also a = and the eigenvalue equation is 



for 



613 = 5„^s-a^-^ 

Pl3 Pl3 



13 



V 7 pi3 



2 Pl3 Pl3 



Fi3 

Pl3 



(72) 



where with n = 1, 2, 5 are symmetric traceless 
3x3 matrices, 5„ is the corresponding polarization five- 
vector for a massive meson with spin 2 (the sum over n 
from 1 to 5 is indicated only by the repeated subscript 
n) . The p-wave wave function P13 and the /-wave wave 
function Fi^ satisfy two coupled equations 



= 



where 



0, 

dp2i f P1ZP2A 
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■-P13" 







P24 



(73) 



G. Triplets with J = 2, or ^7?2 



In this case a — and 



^13 Pl^s'i (pia X 
'-'ii 



Pl3 



P?3 



(81) 



The eigenvalue equation takes the form 



= 



h +k — 

"^sosc \ '^v 

Pis. 



bl3- - I dp24 fpi3P24 W 624 , 
TT 



(82) 



with 



1 + ybis+pL) 



4a" 



P13P24 -^3 



(83) 



2 14, 
Wpp = Jl + a P13P24 ^(3>^2 - -'^o) 

+ (P24 +P13) ^-^1 ' 



p/ 



W/p 



2 

a P13P24 — ^ (3 J2 - Jo) 



45 



V6 



[p24 2Ji-fp?3(5J3-3Ji)] , 

-a P24P13 — (3J2 - Jo) 



(74) 



(75) 



45 



V6 



^ [p^g 2 Ji + (5 J3 - 3 Ji )] , (76) 
16 

W// = (5J3-3Ji)/2 + a2pi3P24^(3J2- Jo) 

- «'(P?3+P24)^(5J3-3Ji), (77) 

and the functions Jo, Ji, J2, and J3, are given in Ap- 
pendix ^ 



F. Singlets with J = 2, or 



In this case ba = and 



, <^ ai3 Plas'^pjs 

<Pl3 — '-'n 9 • 

Pl3 Pl3 

The eigenvalue equation for the function a reads 



(78) 



h +k — 

Pl3 



ai3 / dp24: f P13P24 W 024 

TT 



where 



W = 



1 + y (P24+PI3) 



(79) 



3 1 , , X 
2J2-2^o) . (80) 



V. MASSES AND WAVE FUNCTIONS 

This section describes examples that illustrate to what 
extent the simplest version of the RGPEP approach can 
reproduce masses of the known 66 and cc mesons and 
how the corresponding wave functions may depend on 
the relative momentum of the quarks. 

A. Coupling constant and quark mass 

One potentially valid way to determine the coupling 
constant a\ and quark mass mx in H\ at A = Aq is to 
evolve their values as functions of A using RGPEP from 
the region of large A, say A = Ai, where their values may 
be adjusted to observables that are minimally sensitive 
to the non-perturbative mechanism of binding of quarks 
and gluons. For such observables, the adjustment could 
be made using a perturbative expansion for the 5-matrix 
for quarks and gluons using Hx^ in the femtouniverse ■ 
Although a precisely defined calculation including bound 
states as asymptotic states does not exist yet in the RG- 
PEP approach to QCD, some patterns expected to occur 
in such calculation have already been studied. For ex- 
ample, the RGPEP evolution that starts at Ai must be 
extended down to Aq comparable to the meson mass and 
to reach that far one has to deal with issues of conver- 
gence that require optimization of details of the method 
|53|. It is also known 37] that the differential equations 
of RGPEP that describe the evolution of the operator Hx 
(not the S'-matrix) produce in third-order perturbation 
theory in QCD the coupling constant that evolves with 
A according to the formula 



ao 



CtMz 



1 + [aM^/(6^)] (llA^c - 2n/) In (Aq/M^) 



, (84) 



which matches the well-known formula for the running 
coupling constant in the original Lagrangian calculus for 
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the QCD action |63.l6^. Here, umz is the coupUng con- 
stant in with Ai = Mz, Mz — 92.1 GeV is the mass 
of the Z-boson, Nc = 3 is the number of colors and n f 
is the number of flavors (the theory analyzed here has 
Uf = 1). Thus, as soon as one estimates the value of 
the coupling constant a at one value of A in the boost- 
invariant Hamiltonian approach, such as Ai = Mz, the 
size of a at other values of A is in principle known in the 
entire region, in which the perturbative RGPEP calculus 
for the Hamiltonians can be accurate. 

For example, if one assumes that aMz 0-12 IsJ, Eq. 
(jMjl produces ao 0.326 at Aq - 3.7 GeV (for n/ = 6, 
one would obtain ao ~ 0.21). However, if one uses the 
same formula in a strict expansion in powers of to 
third order only, one obtains a result that is about 30% 
smaller than 0.326. The reason is that Eq. H84I) pre- 
dicts an increasingly rapid growth of the coupling con- 
stant when A decreases. But the perturbative formula is 
replaced by a non-perturbative one for A near the bound- 
state mass, and this can be studied in detail using models 
[soUsij. By analogy with the models, one may expect a 
finite but rapid transition of a from the increasing to a 
decreasing function of A in the region where the binding 
mechanism dominates dynamics, smoothing the disconti- 
nuity present in the perturbative formula (|84|) for a when 
its denominator passes through zero. Therefore, a low- 
order expansion in powers of aMz useful. But it 
is useful to expand the operator H\g in powers of ao- 
The operator coefficients of this expansion can be found 
assuming that ao is infinitesimally small ,52]. The same 
coefficients can be used for evaluating Hx^ when Aq is 
small and ao is comparable with 1. Finding a precise 
formula for ao in terms of a^i for realistic values of the 
coupling constants may require sophisticated high-order 
RGPEP calculations. If the precise value of ao cannot 
be easily calculated in a low-order perturbation theory, 
one can seek values of ao that may correspond to the 
available bound-state data and then incorporate the re- 
sulting picture in a new perturbation theory around the 
first approximation found that way. At the present stage 
of development, one can only verify if the simplest version 
of the boost-invariant Hamiltonian approach can repro- 
duce known masses of heavy quarkonia when the coupling 
constant ao is allowed to take values on the order of 1. 

Less is known about the quark mass m as a func- 
tion of A, and what values of m^o one should expect in 
H\g. Technically, the mass parameter is specified as the 
perturbative eigenvalue of Hx^ for one-quark states 0. 
Therefore, one can expect that the mass should be close 
to the pole mass jSlj, which is about 10% larger than 
the quark mass in the minimal subtraction scheme for 
bottom quarks. One may expect for bb mesons that 
TO = toao ^ 4.5 GeV to 5 GeV. 

Thus, although in principle the boost-invariant Hamil- 
tonian approach appears able to cover the whole range 
of energy scales accessible experimentally in the case of 
heavy quarkonia, one needs to carry out higher-order cal- 
culations than carried out here in order to correlate high- 



energy perturbative results in the femtouniverse with the 
description of binding of quarks and gluons at the scale 
of b or c quark masses in one and the same scheme. In 
the simplest version of the Hamiltonian approach, one 
can only find out if there exist choices for the parameters 
a and m at Ao on the order of the quark masses that 
can produce spectra of masses of the quarkonia with rea- 
sonable accuracy. Since one expects a ~ 1/3, "reason- 
able" means here that the masses should be reproducible 
with accuracy on the order of 1/3 or perhaps 1/10 of the 
largest splittings between states with neighboring quan- 
tum numbers. The latter are on the order of 500 MeV 
and this means that matching data with accuracy on the 
order of 50 MeV would be quite good in the simplest ver- 
sion if the required a and m for such matching are close 
to the values established from other considerations. 

In principle, masses of only two mesons are sufficient to 
fix the values of a and to as functions of A near Ao. The 
question is which two masses one should use. That the 
choice is not obvious and that there is a need for a good 
choice is a consequence of the fact that in approximate 
calculations all masses are calculated with theoretical er- 
rors that are not known and if one uses two masses that 
are obtained with a large theoretical error then results for 
all other masses will be obtained with large errors. Ex- 
perience with exactly solvable models dictates that the 
most accurate procedure should be to choose the masses 
in the middle of the spectrum of the window Hamilto- 
nian 0,|53|- Let us consider the example of bb mesons. 
The most rational choice is to use masses of two p-wave 
mesons xi(lP) and xi(2P). Their masses lie in the mid- 
dle of the window spectrum. The high-energy boundary 
of the window corresponds to short distance dynamics, 
i.e., the most tightly bound states, having the smallest 
masses. The low-energy boundary corresponds to long 
distance dynamics, i.e., the states with largest masses. 
The mesons xi(lP) and xi(2P) are not very sensitive 
to the short distance dynamics and thus also not very 
sensitive to the unknown term a^R in the potential of 
Eq. H49|) , because quarks in these mesons are pushed out 
from the region of small relative distances by the cen- 
trifugal barrier with 1 = 1. One expects that RGPEP 
of 4th-order will produce terms a^R that correspond in 
position space to functions like (5''(r) or Such func- 

tions are known to occur in effective potentials in stan- 
dard dynamics in atomic calculations in QED when one 
includes effects due to the exchange of two photons, ver- 
tex corrections, and self-interactions order a^ [S^. Thus, 
selecting mesons xi(lP) and xi(2P) that have I = 1, one 
has a chance to avoid theoretical errors due to the cur- 
rent lack of knowledge of the terms a^R in LF QCD. At 
the same time, the masses of mesons xi(lP) and xi(2P) 
are most probably less sensitive to the quark-antiquark 
long-distance dynamics than the masses of states with 
I = 2 or 3. At long inter-quark distances, the harmonic 
oscillator potential is expected to lose accuracy because 
the simplest approximation does not take into account 
effective gluons that may be created when quarks move 
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far away from each other. For example, one effective 
gluon could actively participate in the non-perturbative 
dynamics of states with masses that exceed the middle 
eigenvalues of the window Hamiltonian by more than 1 
GeV, which is an estimate of the magnitude of mass of 
an effective gluon at the scale Aq. The estimate indicates 
that one should probably fit parameters a and m in the 
simplest approximation to meson masses that do not ex- 
ceed the middle masses by more than about 1 GeV, and 
the mesons Xi(lP) and xi(2P) lie in this range. 

TABLE I: Qualitative illustration of results of the sim- 
plest approximate approach to heavy quarkonia in the case 
of masses of bb mesons (in MeV). The second column is ob- 
tained using the quark mass m — 4856.92 MeV and couphng 
constant a = 0.32595 for A = 3697.67 MeV when one demands 
that the masses of mesons xi(lP) and xi(2P) are reproduced 
using an auxiliary interpolation procedure described in the 
text, which is employed only to increase speed of numerical es- 
timates in this illustration and is accurate to a few MeV. The 
corresponding oscillator parameters are w = 182.16 MeV and 
kp = 0.157722. The third column quotes experimental data 
with accuracy to 1 MeV and the fourth column displays the 
difference. The fifth column shows precise numerical results 
obtained from the same dynamical equations for the same 
values of m, a, and A, but without errors introduced by the 
auxiliary interpolation procedure. 



meson 


interpolation 


experiment 


difference 


precise 


T10865 


10725 


10865 


-140 


10729.7 


T10580 


10464 


10580 


-116 


10466.9 


T3S 


10382 


10355 


27 


10385.2 


X22P 


10276 


10269 


7 


10278.5 


Xi2P 


10256 


10256 





10258.0 


Xo2P 


10226 


10232 


-6 


10228.1 


T2S 


10012 


10023 


-11 


10013.8 


X2lP 


9912 


9912 


-1 


9913.3 


XilP 


9893 


9893 





9894.2 


XolP 


9865 


9859 


5 


9865.5 


TIS 


9551 


9460 


91 


9551.8 




9510 


9300 


210 


9510.8 



An example of results one obtains by fitting masses of 
bb mesons Xi(lP) and xi(2P) for Aq = 3697.67 MeV is 
given in Table The coupling constant and mass re- 
quired for obtaining the two masses at this Ao are a — 
0.32595 and m = 4856.92 MeV, in a good qualitative 
agreement with expectations (see Eq. (I84|l and the dis- 
cussion that follows it). The large number of digits in 
these numbers is a numerical effect due to the precision 
of experimental data and does not reflect the accuracy of 
the Hamiltonian approach to QCD in its simplest version 
with only \QQ) sector, which is presumably much worse. 
The value of Aq chosen in this example lies in the middle 
of a small range of size of about 200 MeV in which one 
can vary Aq and numerically reproduce the same known 
values of the two meson masses with accuracy better than 
1 ppm by varying the parameters a and m as functions 
of Aq . Table HI shows that a whole set of masses in the 
middle of the window spectrum is also close to data when 



the two selected masses are. The masses near the edges 
of the window, most sensitive to the theoretical errors of 
the simplest version, deviate from data by more than the 
masses in the middle of the window spectrum do, but the 
magnitude of these deviations is not absurdly large. 

Resuhs in Table |I| were obtained in a sequence of steps 
that need to be explained. The key difhculty is that the 
masses can be determined only numerically, and the in- 
tegrals that determine matrix elements of the window 
Hamiltonians depend simultaneously and in nontrivial 
way on a, m, and A. The complication is caused mainly 
by the form factors /, which eliminate the possibility 
of analytic integrations using the oscillator basis func- 
tions. The numerical evaluation of the matrix elements 
of window Hamiltonians takes time. The time required 
for evaluation of one matrix element on a good laptop is 
on the order of a second, and one needs on the order of 
a thousand matrix elements to obtain accuracy of four 
digits for masses of mesons that result from diagonal- 
ization of the window matrix. One would have to carry 
out very long computations to find suitable a and m for 
any given choice of Aq if one were computing matrix ele- 
ments always a new for every change in the parameters. 
Instead, one can evaluate eigenvalues x in Eqs. to 
(|82|l using parameters that lie on discrete points of a grid 
in the parameter space. The parameter kp is more conve- 
nient than A itself. Then one can interpolate between the 
grid points in order to find approximate eigenvalues for 
a continuum of parameters a and m for a whole range of 
values of Aq in the region covered by the grid. Such inter- 
polation produces quickly results of precision better than 
10^'^. The interpolating functions allow one to identify 
the values of a and m that reproduce the same masses 
of xi(lP) and xi(2P) for different values of Aq very effi- 
ciently even though the eigenvalues are less precise than 
to 1 MeV. Results of the interpolation are given in Table 
U in the second column, marked "interpolation." Pre- 
cisely evaluated masses for the parameters selected using 
the interpolation are given in the fifth column in TableQ] 
The quoted values were found stable against a) increase 
of the number of basis states above about 40 (about 20 
if only one orbital angular momentum wave function is 
present), b) changes in the oscillator basis functions of 
Appendix G due to variation in the oscillator frequency 
within about one order of magnitude, and c) changes in 
the algorithm for evaluation of matrix elements of the 
window Hamiltonians (two different integration routines 
produced the same results). Subsequent discussion con- 
cerns results that satisfied the same convergence tests. 

The example given in Table ^ shows that even in its 
simplest version the Hamiltonian approach can lead to 
phenomenologically reasonable results for the masses of 
bb mesons. But the example does not provide information 
about how large is the range of parameters for which the 
simplest version of the Hamiltonian approach can match 
the masses of known heavy quarkonia. This issue is taken 
up using examples in the remaining parts of this section. 
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B. 



Masses of bb mesons 



In order to obtain qualitative information about the 
distance between the simplest version of LF QCD and 
data for masses of bb mesons, one can consider two dif- 
ferent fits of a and m. 

One fit is focused on the middle of the mass spectrum 
of known mesons. Instead of fixing two selected meson 
masses and checking how others are reproduced as it was 
done in the previous subsection, one finds a minimum of 
deviation of the computed masses from data as function 
of a and m assuming different Aq and this is done for 7 
masses in the middle of the experimentally known spec- 
trum. Results of this fit are denoted in Table ^ as "fit 
to middle." They are shown graphically in Fig. ^ Vari- 
ation of the obtained spectrum with Aq when one keeps 
any two of the 7 masses fixed, or rather the degree of 
absence of such variation, is not further studied in this 
and the next subsection. 



TABLE II: Masses of bb mesons (in MeV). The third col- 
umn results from the fit of the coupling constant a and quark 
mass m at the indicated optimal value of A to 7 middle 
masses of known bb mesons, i.e., masses of xo(lP)5 Xi(lP)i 
X2(1P), T(2S), xo(2P), Xi(2P), X2(2P), and this fit implies 
the oscillator parameters co = 184.62 MeV and kp = 0.26667. 
The fourth column results from the fit to masses of all 12 
mesons r;(lS), T(1S), xo(lP), Xi(lP), X2(1P), T(2S), xo(2P), 
Xi(2P), X2(2P), T^Di (estimated as similar to 732), T(3S), T 
10580(84), and the corresponding oscillator parameters are lj 
— 147.11 MeV and kp = 0.016667. Question marks regarding 



D states are explained in the text. 



A [MeV] 
m [MeV] 
a 

T10580 

T ^Di 

T3S 

X2 2P 

Xi2P 

Xo2P 

T 

T ^L»2 

T ^Di 

T2S 

X2lP 

XilP 

XolP 

TIS 

?7fclS 



experiment 



10580 ± 3.5 

10355.2 ± 0.5 
10268.5 ± 0.72 
10255.5 ± 0.72 
10232.5 ± 0.9 
10161.1? ± 2.2 
10161.1? ± 2.2 



10023.3 ± 0.31 
9912.21 ± 0.57 
9892.78 ± 0.57 
9859.44 ± 0.73 
9460.3 ± 0.26 
9300.6 ± 10 



fit to middle 



3779.8 
4835.9 
0.28839 
10734 
10461 
10389 
10273 
10256 
10231 



10115 
10018 
9907 
9892 
9869 
9574 
9542 



fit to all 



3252.3 
4979.7 
0.50738 
10629 
10446 
10329 
10272 
10241 
10194 
10172 
10169 
10154 
9991 
9943 



9849 
9448 
9359 



The other fit is focused on checking how many of the 
experimentally known meson masses can be explained in 
the simplest version of LF QCD, and how accurately. 
The second fit includes masses of all 12 well-established 
mesons |3lj with decay widths significantly smaller than 
100 MeV. A decay width comparable with 100 MeV is 
considered an indicator of relevance of processes that 
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FIG. 1: Illustration of masses in the 3rd column in Table IITI 
The left thick bars in each of the three columns indicate data 
and right thick bars results of computation. The theory mass 
10461 MeV corresponds to a state dominated by the d-wave, 
apparently not easy to identify experimentally [sj. 



the simplest approximate version of LF QCD cannot de- 
scribe. The second fit is denoted in Table ^ as "fit to 
all." The results are shown graphically in Fig. |21 

A separate comment is required concerning the D 
states (d-waves) in Table One such state is known 
experimentally Isij . most probably having total J=2. 
This state is included in the table in order to illustrate 
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FIG. 2: Illustration of masses in the 4th column in Table HTl 
The left thick bars in each of the three columns indicate data 
and right thick bars results of computation. The theory Te- 
state mass changes by only a fe-w MeV, to 10446 MeV, when 
one changes "fit to middle" to "fit to all." 
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what happens in the simplest version of LF QCD con- 
cerning d-wave mesons. All D states, whether one con- 
siders singlets or triplets, or J = 1, are expected in other 
approaches to have similar masses |65l l6^ . The same 
happens in LF QCD. In Table |Il| the state T^Dl is 
known only in theory, and masses of both states T^D2 
and T^Dl are near the one experimentally know mass of 
10161 MeV for the I?-meson that presumably has J = 2. 
The excited triplet state T^i?l is not known experimen- 
tally but comes out of the calculation as dominated by its 
d-wave component. In Fig. ^ for mesons, the known 
mass of 10161 MeV is shown in comparison with the the- 
oretical mass of the T state dominated by d-wave, and, 
in Fig. [3 it is shown in comparison with all three theo- 
retical results for D states with J = 2 and J = 1. 

The range of coupling constants ao and quark masses 
mo for which the Hamiltonians ^qqaq '^^^ match the 
data for bb mesons, as illustrated in Tables III and UTI and 
in Figs. ^ and 121 is summarized using these examples in 
Table jo] (subscript is omitted). Note that the fit to 



TABLE III: Parameters in H\g for bb mesons. 



parameter 


Xi 


7 middle 


all 12 


A [MeV] 
m [MeV] 
a 


3697.67 
4856.92 
0.32595 


3779.8 
4835.9 
0.28839 


3252.3 
4979.7 
0.50738 



all 12 meson masses points to the much larger coupling 
constant and quark mass at considerably smaller Aq than 
in the two similar cases with fits to 7 or only 2 middle 
mesons. This feature most probably emerges because the 
term a^R in the radial factor 1 + a^R in Eq. H49|) is set 
to in the simplest version of the approach. Splittings 
between s-wave mesons, including the rji, and T family, 
are sensitive to the short-distance dynamics that depends 
on the term a^R. Calculation of the term a^R requires 
full 4th-order RGPEP analysis. In the absence of a^R, 
one can nearly reproduce masses of the s-wave states at 
the price of increasing a and to. However, it is clear that 
the 4th-order calculation of the term a^R must be carried 
out in order to narrow the range of possible values of a 
and TO. 



C. Masses of cc mesons 

Masses of cc mesons can be studied in the simplest 
version of the Hamiltonian approach analogously to the 
case of bb mesons discussed in the previous subsection. 
Table IIVI shows results of two fits: to 3 middle masses 
in the known spectrum, and to all 11 masses of well- 
established mesons with small decay widths. Masses of 
two theoretical D states with J = 2 can in principle be 
compared with the one experimentally known mass of 
3836 MeV for a meson whose J = 2 needs confirmation 
[3l|. Fig. 131 illustrates the results for cc mesons obtained 



from the fit to the three middle masses in the window (the 
third column in Table HVjl . the £)- mesons with J = 2 are 
not illustrated. In Fig. ^ illustrating results from the fit 
to the masses of all 11 cc mesons (the fourth column in 
Table llV|l . all three states with d- waves; the triplet J = 1 
state corresponding to ^3770, the triplet state, and the 
singlet state with J = 2, are indicated. 

Table shows examples of parameters that fit masses 
of cc mesons, in comparison to the examples of param- 
eters from Table IIIII that fit masses of bb mesons. It is 

TABLE IV: Masses of cc mesons (in MeV). The third column 
results from the fit of the coupling constant a and quark mass 
m at the indicated optimal value of A to only 3 middle masses 
of cc mesons, i.e., masses of xo(lP), Xi(lP)i X^i^P)^ s-nd this 
fit implies the oscillator parameters uj = 284.93 MeV and kp 
= 3.0642. The fourth column results from the fit to masses of 
aU 11 mesons 7?(1S), J/V(1S), Xo(lP), Xi(lP), X2(1P), '7(1S), 
ip{2S), i/'3770, '04040, ^/>4159, 04415, and the corresponding 
oscillator parameters are cu = 278.72 MeV and kp — 1.3396. 
Question marks regarding D states are explained in the text. 



meson 


experiment ^L] 


fit to middle 


fit to all 


A [MeV] 




1990.0 


1934.2 


m [MeV] 




1553.3 


1577.4 


a 




0.34335 


0.41443 


^4415 


4415 ± 6 


4505 


4462 


^4159 


4159 ± 20 


4178 


4152 


^4040 


4040 ± 10 


4122 


4083 




3836? ± 13 




3801 




3836? ± 13 




3793 


V3770 


3770 ± 2.4 


3773 


3756 


V2S 


3686.093 ± 0.034 


3698 


3662 


77c2S 


3638 ± 5 


3619 


3557 


X2lP 


3556.26 ± 0.11 


3560 


3551 


XilP 


3510.59 ± 0.1 


3507 


3481 


XolP 


3415.16 ± 0.35 


3412 


3340 


J/i^lS 


3096.916 ± 0.011 


3199 


3156 


77c IS 


2980.4 ± 1.2 


3111 


3024 



plausible that the anomalously large result of a ~ 0.5 for 
bottom quarks merely indicates that a fit to all masses 
can push parameters toward explanation of the large s- 
wave splittings at the expense of theoretical consistency 
of the approach. It is expected that a in bb mesons should 
be smaller than in cc mesons, in correspondence with the 
increase of Ao. Calculation of the term a^R in 4th-order 
RGPEP should clarify the situation considerably. 

At this point, one can observe that the Hamiltonians 
■^QQAo obtained for heavy quarkonia in the simplest ver- 
sion of LF QCD seem to explain the tendency of potential 
models to prefer larger values of the coupling constant 
and quark masses than indicated by results based on per- 
turbative QCD j^. Namely, when one forces a single 
non-relativistic Hamiltonian with a Coulomb potential 
at short distances and some confining potential at large 
distances to fit data for masses of many mesons, instead 
of only the middle ones in the window where a simple po- 
tential model can be justified, the strong-interaction rel- 
ativistic effects at short distances are not well described. 
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Similarly, at large distances between quarks, a simple 




FIG. 3: Illustration of masses in the 3rd column in Table HVI 
The left thick bars in each of the three columns indicate data 
and right thick bars results of computation. 




A thorough discussion of potential models in a related 
two-body Dirac formalism is also available [g^. Since 
the Bethe-Salpeter equation or two-body equations are 
ultimately related in QCD to an entire set of the Dyson- 
Schwinger equations [Tfl, one should observe that the 
whole set corresponds in the Hamiltonian approach to 
the eigenvalue problem in which all effective-particle Fock 
sectors are explicitly included. 



TABLE V: Examples of parameters in H\g that fit masses of 
cc mesons, compared with the examples of parameters that 
fit masses of bb mesons. 



parameter 


bb middle 


bb all 


cc middle 


cc all 


Ao [MeV] 
m [MeV] 
a 


3779.8 
4835.9 
0.28839 


3252.3 
4979.7 
0.50738 


1990.0 
1553.3 
0.34335 


1934.2 
1577.4 
0.41443 



FIG. 4: Illustration of masses in the 4th column in Table HVI 
The left thick bars in each of the three columns indicate data 
and right thick bars results of computation. 



A characteristic feature in Table |Vl is that the quark 
mass varies slowly with changes of Aq while the coupling 
constant varies relatively quickly. The width parameter 
A occurs in third power in the oscillator frequency, uj, 
in Eq. H48|) . Therefore, there is a possibility to keep a 
whole set of meson masses approximately constant when 
Ao is changed a little by a considerably larger percentage 
of change in a, while an a priori possible compensating 
change of the quark mass cannot be large because the 
overall scale of masses is dictated by Eq. H54f) . The 
eigenvalues x are negative for the smallest-mass mesons 
and positive for the other mesons. Thus, variation of the 
quark mass is limited by the requirement of preserving 
the relativistic structure of the spectrum in Eq. 1)54(1 . 

Finally, one should stress that the adjustment of a 
and m at A = Ao, which includes a choice of Ao, does not 
provide a check on the renormalization group variation 
of the parameters a and m with A beyond the qualita- 
tive statement of agreement with expectations regarding 
the magnitudes of the parameters. In order to study the 
renormalization group structure, one would have to con- 
sider a plausible choice of a and m at certain Ao, evolve 
the values of a and m in RGPEP for H\ to other values of 
A, and solve the eigenvalue problems for different values 
of A. One would need to include the 4th-order RGPEP 
to begin with and also perform non-perturbative compu- 
tations of the spectra of window Hamiltonians with more 
than one Fock sector built from effective particles. 



potential model cannot reproduce effects due to interac- 
tions that involve gluons "in the air." The parameters of 
potential models have to increase artificially in order to 
keep reproducing the smallest and largest masses of know 
mesons. One should mention that quark models based on 
the Bethe-Salpeter equation with a kernel that resembles 
a harmonic oscillator potential at intermediate distances 
have been introduced a long time ago [b^] and extensive 
studies of the hadronic spectrum have been made using 
models that successfully incorporate such kernels |68| . 



D. Wave functions 

An important aspect of the LF Hamiltonian dynam- 
ics is that it provides wave functions of bound states. 
The discussion that follows is limited to an illustration 
on examples how the Hamiltonian approach works using 
wave functions of mesons J/ip and the ground state of 
T. Masses of these mesons are not described particu- 
larly well in the approximate approach. But their wave 
functions are sufficient to display the main features. The 
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states corresponding to J /tj} and T contain s-wave and 
d-wave wave functions that are invariant under boosts. 
Section III El explains how these states are constructed 
using the wave functions. 

S/k[10-%eV-3/2] s-waves 

1 

bbm bba ccm cca 




1000 2000 3000 4000 5000 
-0.2t k[MeV] 

FIG. 5: Four examples of s-wave wave functions in bb and cc 
mesons. Two curves labeled bb correspond to the ground state 
of T and the two curves labeled cc correspond to J /lii. The 
curves are presented in two versions that correspond to the 
two sets of parameters a and m shown in Tables HTl and HVI 
One set was adjusted to masses in the middle of the known 
spectrum, and the other one to masses of all known mesons 
with small widths. The label "m" refers to "middle" and the 
label "a" to "all." 



D/k[10-6MeV-3/2] d-waves 
1 

Q g [ bbm bba ccm cca 

. 6 
0.4 
0.2 



0.2 




1000 2000 3000 4000 5000 

k [MeV] 



FIG. 6: Four examples of d-wave wave functions in the 
ground state of T and in J/^}). The curves are labeled in 
the same way as in Fig. |S| 

The J /ip and T wave functions have the structure indi- 
cated in Eq. The s-wave wave function is denoted 
by S, and the d-wave wave function by D. Both are 
functions of the relative momentum k of the two effec- 
tive quarks only through its length, k = This is 
a consequence of that the wave functions S and D de- 
pend on the invariant mass squared of the two quarks, 
_ 4(^2 _|_ ;j2-j 'Yhc relative momentum k shares 
many properties with the momentum-space variable typ- 
ically introduced in non-relativistic potential models, but 



one has to remember that the variable k appears in QCD 
according to the rules of LF dynamics. 

The wave functions S and D are shown in Figs. [Sjand 
El in four versions, two for the ground state of T and 
two for J /ip. Two versions per meson are obtained using 
the two choices of parameters a and m that are given in 
Tables im and II VI Numerical values of the wave functions 
can be read from the tables given in Appendix [H] 

It is visible that the c?-wave component is much larger 
in size in the charm case than in the bottom case, al- 
though the s-wave components are similar in both cases 
at small relative momenta. This result can be attributed 
to much more relativistic relative motion of quarks in 
J/tp than in T. Relativistic motion leads to enhance- 
ment of spin-dependent interactions that mix the d-wave 
component with the s-wave component. 

A comment is in order regarding decay widths of the 
mesons. In the leading approximation, the leptonic decay 
amplitudes are proportional to the s-wave wave functions 
at the origin in position space (integrals of the wave func- 
tions in momentum space). There is little doubt that the 
leptonic decay widths in the Hamiltonian approach to 
QCD will be qualitatively similar to the widths obtain- 
able in potential models. On the other hand, inclusion of 
the term a^R in the effective potential and correspond- 
ing self-interactions in HqQ)^, in a calculation similar to 
the simplest version discussed here, and further evalua- 
tion of the gluonic components in the effective-particle 
Fock-space basis in the eigenvalue problem for Hx, may 
provide quantitative insight into effects not accessible in 
potential models. These effects may influence the lep- 
tonic decay rates and need to be evaluated in order to 
precisely compare the theory with data including the de- 
cay rates. The point is that such studies appear feasible 
using the boost-invariant Hamiltonian approach to QCD. 



VI. 



CONCLUSION 



The simplest approximate computation of masses of 
heavy quarkonia in QCD with one flavor of quarks, still 
using the gluon mass-gap ansatz to finesse a simple 
picture in the \QxQ\) sector, suggests that the boost- 
invariant Hamiltonian approach offers a feasible path to 
extended studies of quark and gluon dynamics in the 
light-front Fock space. One can use creation and an- 
nihilation operators for effective particles in order to ex- 
plicitly construct the states of heavy quarkonia. The ap- 
proach produces an approximate constituent picture that 
is relativistic and usable for description of fast-moving 
mesons. Masses of the mesons are reproduced in the sim- 
plest version reasonably well for reasonable values of the 
coupling constant, a, and quark mass, to, using a small 
set of basis states in the effective eigenvalue problem. 
In the simplest version, the short-distance high-energy 
effects and large-distance gluon dynamics are not fully 
described. Therefore, it is not surprising that the sim- 
plest approach can be used to reproduce only a small set 
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of meson masses that lie in the middle of the spectra of 
small effective Hamiltonians called "windows," see Sec- 
tion ^O] But there is also no immediate reason found to 
question that the LF Hamiltonian approach to QCD may 
provide an interesting alternative to other approaches. 

The Hamiltonian calculus produces the boost-invariant 
wave functions that describe heavy quarkonia in terms of 
their virtual effective-particle components in the LF Fock 
space. In principle, these wave functions not only provide 
a relativistic quantum image of a single hadron, but they 
also can be used in description of decays, production, 
and scattering of the quarkonia using QCD. Although 
the cases discussed here concern only bb and cc systems, 
the extension to the case of unequal masses, such as be 
or cb mesons, requires only that instead of the relative 
momentum variable k used here (see Appendix one 
uses the momentum variable defined in analogous way 
by the relations 

= , (85) 
\J rn^ + k"^ + \Jm2 + k'^ 

= R+ZFT^i+Z!. (86) 

Associated momentum-space techniques to handle two- 
and three-particle systems with different masses in the 
context of studies of the bound-state structure or decay 
are sufficiently advanced in the LF approach to handle 
states that contain quarks and gluons j21j]. It is also 
known that the gluon mass ansatz technique works rea- 
sonably well in the case of gluonium [t^. Thus, it seems 
plausible that the case of different quark masses may be 
treated with explicit inclusion of the quark- antiquark- 
gluon sector. Knowing the corresponding wave functions, 
one can attempt to describe a host of new exclusive or 
semi-exclusive processes that involve heavy quarkonia in 
arbitrary motion. 

The formalism of LF dynamics in quantum field theory 
involves a choice of an axis in space, especially in gauge 
theories such as QCD, where one has to make a choice of 
gauge depending on that axis. Therefore, the rotational 
symmetry of the theory is not explicit in the LF Hamil- 
tonian formalism. Most of the expressions one encoun- 
ters depend on the distinguished axis. It is reassuring 
that the LF Hamiltonian approach to heavy quarkonia 
produces in its simplest version developed here explicit 
expressions for bound-state spectra in which masses are 
exactly arranged in multiplets corresponding to the to- 
tal angular momentum (meson spin) J = 0, J = 1, and 
J = 2, and the wave functions of the corresponding states 
are classifiable as waves s, p, d, and /. Nevertheless, the 
complete expressions for the wave functions contain ad- 
ditional relativistic factors that are entirely outside the 
scope of non-relativistic potential models, see Eqs. (|34|l . 
(EHl, (Eni), and Appendix O 

The most attractive feature of the boost-invariant 
Hamiltonian approach to heavy quarkonia, the one that 



makes it an interesting candidate for a new expansion 
method in solving QCD j3|, is that the renormalization 
group procedure for effective particles can be systemat- 
ically studied order by order in expansion in powers of 
the effective coupling constant a\ with A on the order 
of quark masses. Such expansion may provide a rea- 
sonably converging sequence of approximations if a\ is 
much smaller than 1. This study shows that a ~ 1/3 
is a reasonable candidate to reproduce the masses of bb 
and cc mesons in systematic calculations. Genuine 4th- 
order RGPEP studies will further clarify if this hope is 
realistic. 



On the other hand, careful readers have certainly ob- 
served that the LF Hamiltonian dynamics with a har- 
monic oscillator potential leads to the eigenvalues 
that are proportional to the angular momentum in the 
relative motion of quarks, like in the Regge trajectories. 
This is a phenomenologically desired feature, although 
one cannot trust the oscillator picture over large dis- 
tances. But when one considers highly excited states, 
their masses increasing as dictated by the quadratic po- 
tential, the probability of emission of effective gluons will 
be also increasing. A string of gluons may be formed, 
with new potentials between heavy effective gluons that 
require further investigation. The Hamiltonian approach 
could thus lead to a quantum theory of the gluon string 
and provide another reason for the same Regge-like be- 
havior of the spectrum, with a different slope than im- 
plied by the two-quark approximation and with validity 
extending to much larger distances than the size of a typ- 
ical hadron. In fact, for a firm chain of quantum gluons 
to form a string, each pair of the neighboring gluons must 
be held together stronger than by a linear potential, and 
a quadratic potential satisfies this condition. The pilot 
calculation described here suggests that the oscillator fre- 
quencies are on the order of one inverse fermi, and the 
oscillator potential term grows as the relative distance 
squared in fermis with a coefficient given by the quark 
mass. This means that the oscillator potential is strong 
for the inter-quark distances larger than about a fermi 
and the quantum theory of the gluon strings with a sim- 
ilar potential between gluons may turn out to be useful 
in phenomenology. Thus, the structure emerging in this 
pilot study of the boost-invariant Hamiltonian approach 
to QCD has a reasonable chance to grow toward a re- 
alistic physical picture supported by a mathematically 
well-defined theory. This is more than another reason to 
undertake the 4th-order studies of the approach. 



Acknowledgement 



This work was supported in part by the Grants No. 
KBN 1 P03B 117 26 and MNiSW BST-975/BW-1640. 



21 



APPENDIX A: TERMS IN Hx 

This appendix lists details of terms in Eq. (|20() for H\. 
The kinetic energy term for effective quarks reads 



T, 



qX 



E 

(TC 



[k] 



b'' b 



it 

''\kcrc~ 



and for certain A = Aq one can set 



(Al) 



'9^.1' 



(A2) 



with M'^ 



/ X + {k +m ) / {I — x) . The integration 



measure [xk] stands for dx (Pk-^ /[16n^x{l — x)]. That the 
effective mass does not depend on the particle motion is 
a unique property of the RGPEP in LF dynamics. 'rs{x) 
denotes the small-x regularization factors rs{x) rs{l — x), 
where rs{x) — x^ 9{x). The gluon kinetic energy term 
reads, 



E 

(TC 



[k] 



k+ 



(A3) 



but an explicit expression for /i| |36l Is^ is not needed in 
this work. 

The emission and absorption term, Y\ — fxYqgX, is 
= 5E123 /[123]r5(xi/3)r5(x2/3)/A(Xf2,m2) 



j23 ^L^Al^AS + j23 C^L^aWas + h.C. 



(A4) 



where ^23 = ~5 tls 9 tiv ^i* ^ J23 = 5 th 9 3^2 ^i* ^ ^ 
denotes the i5-function of three-momentum conservation 
times IGtt'^, i° with a = 1,...,8 denote 3x3 matrices 
of generators of color SU{2>) gauge transformations for 
quarks, e is the gluon polarization four-vector, j!^^ — 
U2l^U3 and ji^^ = V3l'^V2- 

The potential term, V\ = fxVqqx, is 



Vx 



1234 



[1234] (5t?2<43^A(13, 24) 6J4^4&2 , 



where (see Ref. \^ 
T4(13,24) 
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(A5) 
(A6) 



X [0(2)^5(2:5/1)^5(2:5/4) .F2a(1, 253, 4) 
+ 0(-z)f5(a:5/3)f5(a:5/2)^2A(3, 154,2)] , 

and, for example, 

.F2a(1,253,4) _ xi{Ml2 - m^) + ~ m^) 
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m2)2 + {Mj^ - m'^f 



exp 



A4 



- 1 



(A7) 



A^ij = {ki + kjY is the invariant mass of particles i and 
j. Momenta are labeled according to Fig. [7| P+ denotes 
the sum of plus momenta of annihilated quarks. The sum 



6{z) 



ei-z) 



FIG. 7: Momentum labels in potential terms, 
over gluon polarizations. 



nf^k^ + k^n" 



involves momentum fcg ' = £(2:) {ki ' — k2 



(A8) 



with 



£{z) = 9{z) - e{^z), and 2 = (fc+ - A:+)/(fc+ + k+). 
X5 — \z\, and k^ = k^'^/k'^. The instantaneous interac- 
tion between effective quarks, Zx = fxZqqx, is 



Zx = ~9^Y. / [1234]^t?2i43^A(13, 24) 5^ ^43dA4foA2, 
1234 



(A9) 



where 



Za(13,24) = ^3t23tAh{M%,Ml^) 
X [9{z)fs{xz/i)fs{x^/i) 

+ e{-z)fs{x5/3)fs{x5/2)] ■ 



(AlO) 



APPENDIX B: RGPEP SCALING WITH a 

This appendix includes contributions that originate in 
the (7p,y-parts of the sums over gluon polarizations, which 
were not explicitly described in Ref. |3| ■ 

The analysis of scaling with ao for RGPEP factors 
in the window eigenvalue equation for Hx^ with a mass 
ansatz /i^ is based on the similarity between the structure 
of the eigenvalue Eq. H39|l and the Schrodinger equation 
with Coulomb potential in QED (the same kind of the 
leading picture is also found in Yukawa theory {46] ) . The 
subscripts Aq and are often omitted to simplify nota- 
tion. 

When the relative momentum of electron and positron 
in positronium is written as 



k = 



(Bl) 



where is the reduced mass of the fermions, the 
Schrodinger equation for positronium, neglecting spin ef- 
fects, 



2^ 



V'(fc) 



47r a 



d^k' 

(27r)3 {k-k')'- 



■iP{k') = EtP{k)lB2) 
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takes the form 

/ m - 2 



47r 



_dV 

(27r)3 {p-pj 



.(B3) 



The eigenvalue is i? = a; ^ /i , and the ground-state has 
the eigenvalue E — Eq with a; = xq = — 1 and wave 
function 



00 (p) 



(B4) 



(p2 + 1)2 ■ 

with natural n greater 



Higher states have x 
than 1. 

In the QCD Schrodinger equation with Hqq^ , the self- 
interaction terms and the potential kernel contain sim- 
ilar expressions. The self-interaction terms are easy to 
analyze if one knows how to analyze the structure of 
v\ in the eigenvalue Eq. At certain A = Aq, 

with Aq parameterized according to Eq. 1)42(1 . one has 
z;ao(13,24) = z;o(13,24) and 



expressions for (1, 253, 4)/A:^ and J^{3, 154, 2)/k^ that 
contribute to V using identities 



M 



m 



253 



Xi 



AI53 - TO^ 

X4 



-M24, 



(B12) 



54 



3^3 



X2 



24 • 

(B13) 



One starts with expressions, see Eq. (IA7|) . 
J^2(l,253,4) _ 1 2;i(Mi2-n^^)+2:4(Mi3-"i2) 



fc5^ 



■^2(3,154,2) 
'^5 



where 



'53 



X (// - 1) , (B14) 
_ 1 X2 (TW^i - TO^) + (A1^4 - m^) 



(Mil - m2)' + (A^i4 - m2)' 



X (//-I) , 



(B15) 



Wo 
B 



i 12^43 
p+2 



where 



1 

N 
1 d 

--.2 I 



/o(13,24)-^ +i 1^0(13,24) 



/o(13,24) 



V 



1 



wo(13,24) 



^/o(13,24)Z, 



+ 



+ 



e{z) f5(x5/i)f5(a;5/4) /52/! 
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|z|(m2-Ali3)/a;4-M'(2,5,3) 
6>(-z) fj(x5/3)rj(a;5/2) /54/51 
|z|(m2-X2^)/x2-/i2(l,5,4) ' 



(B5) 
(B6) 

(B7) 



V = 0(z)f5(x5/i)r5(x5/4).F2A(l,253,4) 

+ 0(-z)f5(a:5/3)^5(a^5/2)^2A(3,154,2), (B8) 

Z = 0(z)f5(a;5/i)fa(x5/4) 

+ ^(-z)f5(a;5/3)fa(x5/2) , (B9) 

Wo _ 6l(z) r^(x5/i)fa(a;5/4) 752/53 

\z\ - |z|(m2-Mi2)/xi-M2(2,5,3) 

6'(-2:) r^(a;5/3)^.5(a^5/2) /54/51 

^ |z|(m2-Xi4)/x3-/i2(l,5,4) 



Lz ^ ^ V 2x1 



-^53 



2X4 



V 2X2 



2X3 



(BIO) 



(Bll) 



and /y denotes f\o[m'^, {h + %)^]- 

In order to describe the structure of vq for relative 
quark momenta comparable with the strong Bohr mo- 
mentum, introduced in Eq. H41|) , it is convenient to write 



and 



Mil - 



M 
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Mi, - 

MI4 - 



= Vi 
- ^ 

= 2?3, 

= ^V: 
X4 



X]_ 

1^1 
1^1 



'*'13 



The definitions include 



Xi 

z 

— I 

X3 



Xi - X2 , 



(B16) 
(B17) 

(B18) 
(B19) 



(B20) 
(B21) 

(B22) 



and it is helpful to use three-dimensional CMF relative 
momentum variables fcis and ^4, and q — fci3 — ^24- So, 
for ij = 13 and 24, 



7^3 



and 



'''13 



«,24 



M 



13 



24 



(B23) 



(B24) 



The first step is to establish that the potential does not 
generate any small-x singularities in its fully relativistic 
form 01 . The next step is to analyze scaling with a. 
The key to scaling with a for given quark mass m is the 
substitution 



fcs Pi J , 



(B25) 
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where fc^ is the strong Bohr momentum of Eq. (|41|l . 
The dimensionless variables pij, with ij = 13 or 24, are 
typically on the order of 1 in both the purely Coulombic 
case of QED and in the QCD case that includes the har- 
monic oscillator potential studied here. A dimensionless 
momentum transfer p is defined by 



ks P, 



(B26) 



so that p = pi3 — P24- Factors // limit \p\ to values or- 
der (4a/3)^'^Ap, and the additional damping due to e > 
provides a possibility to formally separate the dominant 
terms in the limit a —^ because the Coulomb eigenvalue 
problem is dominated by the dimensionless momenta pij 
on the order of 1. The outer- most factor /g in the poten- 
tial terms limits changes of momenta from above by 
Q,€-o.5 g^j^^ -(-j^jg becomes irrelevant for very small a, 
leaving the Coulomb interaction and the harmonic oscil- 
lator term that provide the leading approximation. 
Observe that 



■^2(1,253,4) ^ (Mi53-C253)k|-^ 

■^2(3,154,2) ^ iMl,i-Ci5i)\z\-' 
k+{f.f-l) Mt,^-2Ml,^Ci5i + Di54 



where 



(^2 5 3 — 



■D253 — 





\- xlMli 






xl- 




xlM\, - 


h xlM^i 






xl- 


Vx\ 


xjMj,- 


V x\M\^ 






4- 


Vxl 




V x\M.\^ 






4- 


Vxl 



C154 — 
-D154 = 

Using Eq. (|B25(I . and introducing two three- vectors 



,(B27) 
,(B28) 

(B29) 
(B30) 
(B31) 
(B32) 



I* fci3 + /S24 ^1 ^ 

m 



one obtains 



■^253 ~ '^253 — -^154 ~ ^"154 



(B33) 
(B34) 

(B35) 



0(a5) 



^253 - <-^253 — -'^154 " 0154 H — 



(B36) 
(B37) 
(B38) 



Thus, 



.^2(1,253,4) ^ .^2(3,154,2) 

kt{!]-^) ~ k+{ff-l) 



+ Oia) (B39) 



q 



\ + ^^^i^^^-fO(a). (B40) 



The effective-gluon exchange term wq in Eq. (jB10|l . is 

2e{z) 7=5(0:5/ 1)^5(3:5/4) 752/53 



Wo 

\z\ 



g2 + ^2(2, 5, 3) + 0(^3+4.) 
2e{-z) fs{x^/^)fs{x^/2) /54/51 



(B41) 



g2+;i2(1^5^4)+0(^3+4.) 

and the intermediate gluon spin sum contributes 

d = q^ + 0{a^). (B42) 



In summary, the factors A and B, defined in Eqs. ljB6|) 
and IB 71). scale as 



A 
B 



—fa 
-fo 



1 

q 



qi 



1 (f '^^ 

// (. q' 



+ c 



e^r - e^r -l- 0{a ) , 



,(B43) 

,(B44) 
(B45) 



and c ~ q;2 because |^ | = 0{a). These scaling results 
are valid even if the mass ansatz /i^ is of the order of a 
instead of 1. 

Scaling analysis of the self-interaction terms begins 
with the RGPEP expression for renormalized effective 
quark mass terms in the eigenvalue Eq. H39|l . 




where 



M 



y 



y(i - y) 



(B46) 

(B47) 
(B48) 



and ji^ is the mass ansatz for the effective gluon accompa- 
nying the quarks i' and j . One can introduce the variable 



q^=p' 



ym, 



(B49) 



and observe that when Ao ~ a'^-^'^'^m, the magnitude of 
q is limited by the RGPEP form factor to the range 
between and a^'^^'^m. Corrections order y cancel out 
or multiply terms order . If /i^ jg of order 1, deviations 
from y = introduce corrections order a^"*"^' and can 
be formally neglected when one keeps only terms order 
and a^. Note that the terms order cancel out 
completely, and even for order a the first correction 
due to y 7^ 1 is order a'^'^^'^. Then, since the integrand 
is symmetric in t = cos0, where qz = qt, and since xiX2 
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differs from 1/4 in order a^, the leading contribution from 
the self-interactions is 



Xi 



X3 



8m-- 



3 167r3 



X // ( ^ - ^ 



(B50) 



APPENDIX C: SPINORS 



The 4x4 matrix wave function in Eq. 1)35(1 . is 
written using spinors m^.s and Vk_s that are obtained by 
applying matrix -B(fc, to) from Eq. H37|l to spinors at rest, 
uq.s and uo,s, which are defined as 



uq^s = v2to 
wo,s = 2s-\/2to 











(CI) 
(C2) 



Xs denotes standard two-component spinors, with upper 
component equal 1/2 + s and lower 1/2 — s, and s means 
— s. Instead of s = ±1/2, u = 2s with values ±1 is often 
used below. The above spinors correspond to fermions 
at rest in the frame of reference in which one carries out 
the calculation for the bound state with momentum com- 
ponents P+ and . The discussion below is simplified 
to the case = since the transverse motion of the 
bound state does not introduce any change in the final 
formulas. 

The wave function "^{kij, Si, Sj) in Eq. H36|l is defined 
using the matrix '^cMFij{kij) and spinors w^:..^. and 
v_j: ^ that describe the fermions in the CMF of the 
quarks i and j. These are obtained by standard matrices 
for boosts along kij instead of the LF boosts. Namely, 



m 



L{k)V2m 

L{k) ct\/2to 
1 



X-kij 






(C3) 
(C4) 

+ I3m\l3. (C5) 



The two component spinors at rest are turned away from 
the z-axis using a 2 x 2 matrix C(^): 



(C6) 



k+ y/2m{Ek + to) 
m {k+ + mf + fcJ 



(C7) 



The matrix ^ introduces the spinor basis in which the 
wave function '^CMFijikij) satisfies a rotationally sym- 
metric eigenvalue equation in the leading approximation. 
The matrix C has been used before by Melosh ,42,] as a 



candidate for description of how the constituent quarks 
are related to current quarks, and by Brisudova and 
Perry 49] in LF eigenvalue problems. Here, the con- 
stituent quarks are dynamically related to current quarks 
using RGPEP and C plays only the kinematical role of 
choosing a basis for spinors. The matrix C is needed be- 
cause there is a change of frame of reference involved in 
expressing LF spinors in the frame of reference in which 
the whole quarkonium has momentum P+ and P^ — 
in terms of the spinors in the CMF of the pair of quarks. 
For example, Wfci.si = ^^s'^ki^^si^ where 



= A^ 



(C8) 



When one uses the slightly rotated basis for the two- 
component spinors in the CMF of fermions, as indicated 
in Eqs. I|C3|) and l|C4|l . and than calculates the spinors in 
the frame of reference where the bound state eigenvalue 
P~ is calculated, one obtains spinors that are used in Eq. 

m- E.g., 

Ui = CijB{hj,m)C^{k,j)L-^{k,j)u^^^^^^. (C9) 
The matrix Q is defined to render 

B(fci3,TO)r'(fcl3)isD(fcl3) = 1, (CIO) 



and Ml = £i3U^^3 ^,^. Similarly, U2 = £24 wj^^^^^, 

-^24 ,S4 ' 



-^13 «_fci3,s3; and VA = £24 w_t„ 



APPENDIX D: BREIT-FERMI TERMS 

In the leading approximation for small a, the potential 
kernel given in Eq. (|40|) is 



Wo 



J ^2 J12J43 



4to^ ^12^43 gpt^ ^^2^43 



,(D1) 



The Breit-Fermi terms in this article originate from the 
first term. Note that the first term contains only one 
form factor and this form factor limits the change of the 
square of the dimensionless momentum pij in the scal- 
ing analysis by a number of the order of a^'^~^ , which is 
much larger than 1 when < e < i, see Appendix IBl In 
contrast, the second term contains two form factors and 
in the scaling analysis these form factors limit the dimen- 
sionless momentum transfer p = Jp ] between quarks by 
a small number on the order of {\pz\/p)X^{'^a/?>Y'^ , see 
Section IIII Dl and Appendix El This difference between 
the form factors implies, that in the first term the domi- 
nant momentum scale is of order 1, originating from the 
Coulomb potential, while in the second term the allowed 
momentum exchange p is in principle infinitesimal as long 
as e > 0. Thus, the second term would be negligible in 
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the scaling analysis if it did not contain the diverging 
factor . This divergence is regulated by the falloff of 
the ansatz as a function of qz when tends to zero, 
as required by the condition that the ansatz does not 
produce a small- a: divergence 

In the leading approximation, the current factors are 
diagonal in spin: ji2i'&ti equals 4to^, and jj^j^Jj always 
equals 4P+^ times ,^xiX2X^X4^. The square root reduces 
to 1/4 since all the xs differ from 1/2 only by terms of 
order a or smaller. In this case, the second term in Eq. 
(|D1|I is the same as the integrand in Eq. ljB50p for self- 
interactions. The self-interaction and the second term in 
Eq. (jDl|l produce together the harmonic oscillator po- 
tential in Eq. H5()|l with the dimensionless spring constant 
given in Eq. ((SSI 4]. 

Beyond the leading order, one has to analyze the factor 
j 12 j 4:3^1- It can be re- written in a matrix notation of Eq. 
(|51|l . The BF terms in the potential kernel V act on the 
2x2 matrix wave function (p which is defined as follows. 
Using results from Appendixinj the sum over quark spins 
in Eq. H39|) can be written as 



S2S4 (^24) 



(D2) 



= {jti-j +m) '^cMFijihj) {h] ^ttt) : (D3) 



13 vt. 



where k = {k^,k), k = {k^,—k), and fc° = vm^ + fc^. 
This 4x4 matrix notation can be replaced by a 2 x 2 
matrix notation using 



where 



(D5) 



is the wave function that appears in Eq. I|5U|) : a and b are 
together four functions of the dimensionless relative mo- 
mentum pij . The fourth root in front of (j) is introduced 
because the measure factor, 



dX2(Pn2i 
X2X4 



2(fiko 



^Jn? + kl^ 



(D6) 



needs to be symmetrized with respect to relative mo- 
menta ^24 and fci3. The resulting potential contains a 

product of ^1 + k^.^/m? and ^J'\~+^^^Jm^ in denomi- 
nator, and the integration measure becomes (fik24 like in 
a non-relativistic Schrodinger equation. Thus, the lead- 
ing contribution of the measure to BF terms through M 
in Eqs. (USjl and jSll) is 



The spin contribution 1 + a^S in Eq. H49() is obtained 
from the sum over quark spins 



S2S4 



4m2 



1 -I- k2A,/rrfi 



(D8) 



using the two component spinors in Eq. ljC6p and the 
wave function in Eq. IjDSp . One multiplies the whole 
eigenvalue equation by xi from the left and by (T3X3 from 
the right and sums up over spins 1 and 3. Then, the 
kinetic energy multiplies only the matrix 4>{pi3), and the 
potential term contains the matrix 



(D9) 



where, using = 7^7^ and notation from Ref. [t^ . 

Si" = [^^^+"'^^^"]/:i3a^/:24 , (Dio) 

CM \-k24,a,E24+m] /• u n -kiaff /t-vi i ^ 



When one neglects terms that vanish faster than in the 
scaling analysis, the matrices £13 and £24 are equivalent 
to 1 and the resulting matrix in the potential takes a fully 
rotationally symmetric form. 



E24+m 

^24'? 



\—k24ff,E24+rn\ 



Ei3+Tn 



■im^{Ek^3 +m){Ek24 + m) 



The result is 4- (4a/3)2S'/16, where 



(D12) 



S = (P?3 + P24)<?^ + Pl'i^ P24CT + P2iS Pl3^ 
+ (T* P2aS4> P24CTCr* + cr* P2A^4' P13S 
-f pi3(7 CT* P24(?cr* +pi3(7Cr' cr* pi3(T. (D13) 

The first term in S is canceled by M from Eq. ljD7|) , and 
after summing over i — 1, 2, 3, one obtains Eq. H55|l . 
Useful identities for Pauli matrices include 



a' baa' = -bS , (D14) 

a'aabScr' = a a b a + 2b(7 aa , (D15) 

a' a abac a a' = ~aa ba ca + 2ba ca aa 

- 2aacdba. (DIG) 

APPENDIX E: ANGULAR INTEGRALS 

The generic form of the integrals over angles in Eq. 
(EOJ is 



A'/ = (^a) {P?3+Pi4) • (D7) 



pi- 



q'^q^ ...q^ 



(El) 
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Using 



Jn. — 



I + — 2pqz ' 



one has 



p2 q2 
2pq 



Jn = 1 Jn-1 



2pq n ' 



and 



r 

jij 
jijk 



27rJo, 
27rg Ji Bp , 

7rq^ [{Jo~J2)s;^ +2J2e;el,] , 
TT ^3 (J, ^ jg) (e^ + + s]^) 



where 



Jo 
Ji 



and 



M \p-q\' 
jJq , 

pq 

l'^ Jo 7 , 

pq 

l'Jo--(l' + \ 
pq \ 3 



= 6'^ - el 



(E2) 
(E3) 

(E4) 

(E5) 
(E6) 
(E7) 

(E8) 

(E9) 
(ElO) 
(Ell) 
(E12) 



(E13) 
(E14) 



APPENDIX F: BASIS FUNCTIONS 

The kinetic energy and harmonic osciUator interaction 
term are of the same form in the eigenvalue equations for 
all mesons, 



2Uosc P ^p ; 



(Fl) 



with the spring tension fcp given in Eq. 153(1 . The Hamil- 
tonian that provides the basis for solving the eigenvalue 
equations is 



2Hb ^ p^ ~ kb Ap 
The eigenfunctions of 2Hi,, 

(f'nlmip) = 4'nl{p)Yim{^p) , 



(F2) 



(F3) 



2/ + 1(1 ™)! pm(cos0)e""*, (F4) 



TABLE VI: Coefficients s„ and d„ in 
of the ground state of T in two cases 
columns third and fourth in Table ITTl 



Eq. IGU in the case 
corresponding to the 
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fit to middle 


fit to middle 

d„ 


fit to all 


fit to all 


1 


0.91359 


0.00497 


0.73034 


0.00532 


2 


0.33678 


0.00462 


0.48078 


0.00627 


3 


0.17577 


0.00403 


0.33108 


0.00632 


4 


0.10675 


0.00344 


0.23735 


0.00595 


5 


0.07010 


0.00289 


0.17449 


0.00539 


6 


0.04807 


0.00239 


0.13013 


0.00473 


/ 


U.Uoo ( o 


n nni op; 
U.UUiyo 


u.uy (to 


U.UU4U0 


8 


0.02405 


0.00156 


0.07356 


0.00340 


9 


0.01723 


0.00123 


0.05533 


0.00280 


10 


0.01237 


0.00096 


0.04150 


0.00226 


11 


0.00887 


0.00074 


0.03101 


0.00180 


12 


0.00635 


0.00056 


0.02308 


0.00141 


13 


0.00453 


0.00042 


0.01710 


0.00110 


14 


0.00322 


0.00031 


0.01262 


0.00084 


15 


0.00228 


0.00023 


0.00928 


0.00064 



contain the radial wave functions 4'ni (p) that satisfy 
l{l + l)kb 



-dpP dp 



Xb 



<PM = 0.(F5) 



In terms of the scaled variable q = Plk\^^ and the 
eigenvalue Xb — yVkb, the substitution (f){p) = x{q)/<1 
produces 

- X" + X{q) +q'x ^ y X{q) ■ (F6) 

q^ 

Eigensolutions normalized to J dqq'^ \ip{q)\'^ — 1 are {L 



TABLE VII: Coefficients s„ and d„ in Eq. (IGU in the case 
of J/il^ in two cases corresponding to the columns third and 
fourth in Table HVl 





fit to middle 


fit to middle 


fit to all 


fit to all 


n 




d„ 






1 


0.94793 


0.01593 


0.91858 


0.01879 


2 


0.27558 


0.01302 


0.33242 


0.01631 


3 


0.13007 


0.00993 


0.17069 


0.01285 


4 


0.07212 


0.00723 


0.09910 


0.00957 


5 


0.04241 


0.00504 


0.06014 


0.00680 


6 


0.02540 


0.00338 


0.03695 


0.00463 


7 


0.01524 


0.00218 


0.02268 


0.00305 


8 


0.00910 


0.00137 


0.01386 


0.00196 


9 


0.00541 


0.00085 


0.00843 


0.00124 


10 


0.00321 


0.00052 


0.00511 


0.00078 


11 


0.00190 


0.00032 


0.00310 


0.00049 


12 


0.00113 


0.00019 


0.00188 


0.00030 


13 


0.00067 


0.00012 


0.00114 


0.00019 


14 


0.00040 


0.00007 


0.00069 


0.00012 


15 


0.00024 


0.00004 


0.00042 


0.00007 



47r {l + my. 



denotes generalized Laguerre polynomials and P{n, k) 
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Pochhammer symbols) 

y = 4n + 2; + 3, 

Xni{q) = (-1) 



2 n! 



r(n + ; + 3/2) 



r(A + n+ 1) 



Tin 



P{~n,k) = nf-io (-n + m) 



1) ^ 



fe=0 



P{—n, k) 
T{\ + k + l)k\ 



(F7) 



(F8) 
-(F9) 



(FIO) 



The oscillator eigenvalues and corresponding radial basis 
functions in momentum space, normalized to 1, are 



Figs. 13 and El Tables IVll and IVIll contain first 15 coeffi- 
cients Sn and dn in the expansion of functions >5'(fc) and 
D{k) in the basis introduced in Appendix IfI with fcf, — 
kp. The first 15 coefficients arc sufficient to see how fast 
the expansion of eigenstates in the used basis converges. 
The actual calculation of these coefficients involved 30 
s-wave and 30 d-wave basis states. The precision with 
which the coefficients are produced in Tables VI and VII 
is limited to the digits that were stable when the num- 
ber of basis states used in the calculation was increased 
above about 20 per wave, fc^ is the Bohr momentum. 



Xb = {An + 21 + 3) ^/h , (Fll) 
Mp) = Xniip/kl^*)-^. (F12) 



S{k) 
D(k) 



{k/ks) 



(Gl) 



APPENDIX G: DETAILS OF THE WAVE 
FUNCTIONS 



This appendix provides numerical data concerning the 
wave functions S{k)/k and D{k)/k that are shown in 



The wave functions are normalized to / dk{S^+D^ 
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